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ABSTRACT/SUMMARY 


Acoustic-mode  combustion  instability  has  long  plagued  the  solid- 
propellant  industry,  and  the  increasingly  frequent  requirement  for  "reduced- 
smoke"  propellants,  with  concommitant  removal  of  metals  from  the  propellant 
formulations  (and,  thus,  removal  of  metal  oxide  particulate  products,  which 
have  a  major  role  in  damping  of  acoustic  oscillations)  is  expected  to  exacer¬ 
bate  this  problem.  One  strategy  for  alleviating  the  problem  involves  identi¬ 
fication  and  utilization  of  approaches  to  decreasing  a  major  source  of 
acoustic  energy,  namely,  the  transient  burning  rate  response  of  the  solid 
propellant  to  pressure  and/or  crossflow  velocity  oscillations.  Previous 
preliminary  modelinn  studies  have  indicated  that  it  might  be  possible  to 
decrease  the  pressure-coupled  response  functions  of  composite  propellants  by 
suitable  modification  of  the  relative  activation  energies  of  the  fuel  and 
oxidizer  ablation  processes. 

However,  the  modeling  utilized  in  arrival  at  this  preliminary  con¬ 
clusion  was  very  oversimplified;  accordingly,  the  purpose  of  the  current  study 
was  to  examine  this  concept  further  using  a  much  more  sophisticated  composite 
propellant  combustion  model  (steady-  and  unsteady-state).  During  the  course 
of  this  work,  such  a  model  was  developed  and  exercised  parametrically  to 
define  potential  effects  of  varying  either  oxidizer  or  fuel  ablation  energy  on 
the  all-important  real  part  of  the  pressure-coupled  response  function. 
Results  of  this  study  indicate  that  additional  factors  not  included  in  the 
previous  simplified  modeling  greatly  diminish  the  hoped-for  beneficial  e'^'fects 
of  increasing  the  activation  energy  of  the  fuel  ablation  process  at  a  constant 
value  of  oxidizer  ablation  activation  energy.  With  respect  to  variation  of 
the  activation  energy  of  oxidizer  ablation  at  a  fixed  value  of  the  fuel  abla¬ 
tion  activation  energy,  the  previously-calculated  beneficial  effects  of 
reducing  oxidizer  ablation  activation  energy  are  strongly  reversed;  in  addi¬ 
tion,  increasing  the  activation  energy  from  the  baseline  value  associated  with 
unmodified  ammonium  perchlorate  has  only  a  min''scule  beneficial  effect. 
Accordingly,  the  initial  conclusion  of  this  study  is  that  the  originally 
postulated  approach  to  decreasing  pressure-coupled  response  functions  of 
composite  propellants  via  suitable  modification  of  the  oxidizer  and/or  fuel 
ablation  activation  energies  actually  holds  little  promise. 
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2.0 


INTRODUCTION/BACKGROUNO 


Over  the  years,  the  appearance  of  combustion  instabilities  in  solid 
rocket  motors,  in  developmental  stages,  and  even  in  operational  systems  (as 
the  propellants  age,  with  changes  in  burning  characteristics)  has  cost 
hundreds  of  millions  of  dollars  in  developmental  and  refurbishment  costs  as 
well  as  leading  to  performance  compromises  associates  with  "fixes"  to  these 
instabilities.  (Such  acoustic  instabilities,  occurring  over  a  wide  range  of 
frequencies,  are  generally  intolerable  in  that  they  can  lead  to  unacceptable 
interactions  with  control  systems,  structure,  and/or  payload  and,  in  some 
cases,  to  catastrophic  motor  failure.)  Before  approximately  1960,  insta¬ 
bilities  occurred  quite  frequently  in  solid  motors,  particularly  tactical 
motors  where  natural  resonant  frequencies  are  typically  in  the  1  Khz  and 
higher  regime.  A  fortuitous  "fix"  to  many  (though  not  all)  of  these  insta¬ 
bilities  was  found  to  be  addition  of  large  amounts  of  metal  additives  to  the 
propellant  formulations,  with  the  resulting  condensed-phase  products  producing 
large  amounts  of  damping  of  oscillations  due  to  velocity  lags  dissipating 
large  amounts  of  acoustic  energy.  (In  fact,  large  amounts  of  metal  were  first 
added  to  propellants  just  for  this  effect,  before  it  was  found  that  such  metal 
addition  also  provided  significant  performance  gains.) 

However,  increasingly  stringent  mission  requirements  (with  increased 
energy  density)  have  led  to  increased  incidence  of  instabilities  in  modern 
rocket  motors,  even  with  large  amounts  of  metal  additive.  In  addition  (and 
probably  most  important),  there  is  a  major  recent  shift  to  a  requirement  for 
"reduced  smoke"  propellants  for  tactical  missiles  to  improve  survivability  (by 
reducing  detectability)  of  the  missiles  and  their  launch  platforms,  e.g., 
aircraft.  This  requirement  precludes  inclusion  of  more  than  very  minor 
amounts  of  metal  particles  in  propellant  formulations,  leading  to  a 
considerable  increase  in  the  occurrence  of  instability  problems  due  to  removal 
of  a  major  acoustic  energy  sink  (damping  term)  associated  with  the  condensed- 
phase  products.  Accordingly,  it  is  important  that  an  alternative  approach  to 
alleviating  instability  problems,  namely  reduction  of  driving  (acoustic  energy 
source)  terms  be  addressed.  One  major  source  term  (generally  the  most  impor¬ 
tant  source)  is  associated  with  interaction  of  an  acoustic  wave  with  the  solid 
propellant  combustion  processes.  Both  pressure  and  flow  oscillations 
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associated  with  an  acoustic  wave  can  cause  solid  propellant  burning  rate 
(energy  release  rate)  oscillations  which,  in  turn,  can  feed  energy  into  the 
acoustic  wave. 

The  relationship  between  burning  rate  oscillations  and  pressure 
oscillations  is  referred  to  as  the  pressure-coupled  response  function  (a 
complex  variable)  of  the  propellant,  while  that  between  burning  rate  oscilla¬ 
tions  and  velocity  oscillations  is  referred  to  as  the  velocity-coupled 
response  function  (although  the  latter  is  generally  agreed  not  to  be  an 
intrinsic  propellant  property,  depending  strongly  on  the  structure  of  the 
local  mean  flowfield).  Attention  here  will  be  restricted,  at  least  for  the 
time  being,  to  the  pressure-coupled  response  function;  it  is  easily  shown  that 
the  real  part  of  this  complex  response  function  is  the  important  quantity  as 
regards  driving  of  acoustic  waves,  with  relationships  between  the  growth 
constant  for  oscillations  and  the  pressure-coupled  response  function  being 
presented  and  discussed  in  numerous  references,  most  notably  CPIA 
Publication  290,  "Combustion  Instability  in  Solid  Rocket  Motors. For 
typical  scenarios,  it  is  found  that  the  pressure-coupled  response  is  a  major 
contributor  (driver)  in  the  acoustic  energy  balance. 

Accordingly,  development  of  means  of  reducing  the  amplitudes  of  the 
pressure-coupled  response  functions  of  a  wide  range  of  solid  propellants, 
leading  to  significant  reduction  of  the  driving  terms  in  acoustic  energy 
balances  for  motors,  should  lead  to  major  reduction  or  elimination  of 
instability  problems  in  motors  utilizing  such  propellants.  Thus,  it  is 
apparent  that  development  of  a  systematic  approach  (principle)  to  reduction  of 
the  magnitude  of  composite  propellant  pressure-coupled  response  functions 
would  be  of  tremendous  potential  importance  to  the  solid  rocket  field. 

For  many  years,  analysts  have  attempted  to  estimate  pressure-coupled 
response  functions  for  composite  propellants  using  theoretical  approaches  (in 
general,  equivalent  to  the  well-known  ZN  approach) which  are  strictly 
applicable  (at  best)  to  homogeneous  propellants.  As  shown  by  King^^"^)  in  the 
early  1980's  (utilizing  very  simplified  composite  propellant  combustion  models 
to  demonstrate  the  point  clearly),  composite  propellants  provide  an  additional 
potential  destabilizing  mechanism  beyond  those  treated  by  the  ZN-type 
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approaches.  This  mechanism  involves  production  of  oscillations  in  the 
oxidize'^/fuel  ratio  of  "packets"  of  propellant  gases  leaving  the  surface  of 
the  propellant  (and,  thus,  oscillations  in  the  gas-phase  heat  release)  accom¬ 
panying  oscillations  in  surface  temperature  due  to  different  activation 
energies  for  the  fuel  and  oxidizer  ablation  processes.  A  brief  outline  of  how 
this  comes  about  follows. 

The  oxidizer/fuel  ratio  of  each  "packet"  of  gas  leaving  the  com¬ 
posite  propellant  surface  may  be  calculated  as: 


OF 


A^exp(-Eox/RT3)  S, 


ox 


packet  A2exp(-E^Jg^/RT^) 


(1) 


where  and  are  the  nondimens ional  surface  areas  associated  with  oxidizer 
and  fuel  respectively.  At  steady-state,  overall  continuity  forces  these 
surface  areas  to  adjust  such  that  OF  is  equal  to  the  overall  propellant 
oxidizer/fuel  ratio,  independent  of  the  mean  surface  temperature,  T^.  Thus, 


prop 


A2exp(-Efugi 

AiexpV'E^/'RT-sT" 


(2) 


Now  if  under  oscillatory  conditions  the  frequency  is  sufficiently  high  that 
Sq^  and  cannot  adjust  significantly  away  from  their  steady-state  values,  it 
may  easily  be  shown  through  combining  of  Equations  1  and  2  (linearizing 
around  T^)  that: 

“fpacket  ■  %rop'’'P[<Eo.-  Ef„el)IV  (3) 

(It  is,  of  course,  recognized  that  in  actuality  there  will  be  oscillations  in 
surface  structure  at  least  partially  compensating  for  the  differential 
temperature  sensitivities  of  fuel  and  oxidizer  ablation,  with  the  degree  of 
compensation  increasing  with  decreasing  oscillation  frequency  until  total 
compensation  occurs  as  the  frequency  approaches  zero;  this  is  treated  in  the 
current  modeling  activity  as  discussed  in  later  sections  of  this  report.)  For 
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a  typical  AP/HTPB  composite  propellant,  the  available  chemical  enthalpy  is 
related  to  the  oxidizer/fuel  ratio  approximately  (again  linearizing)  by: 

"packet  ■  3.87)  (4) 

leading  to: 

"packet  ■  S'"’  ♦  180(0Fp^„pexp[(E„^-Ef„^,)(T^-T^)/RT^2,  .  33^^  ,5, 

Based  on  the  current  literature,  it  appears  that  the  ablation  of 
ammonium  perchlorate  has  a  somewhat  higher  activation  energy  than  the  ablation 
of  typical  HTPB  binders  (by  approximately  5  kcal/mole).  Thus,  upward  pertur¬ 
bations  in  surface  temperature  will  be  accompanied  by  similar  perturbations  in 
oxidizer/fuel  ratio  of  "packets"  of  gas  leaving  the  propellant  surface;  since 
AP/HTPB  formulations  are  generally  fuel-rich,  this  will  in  turn  lead  to 
increases  in  gas-phase  heat  release  (following  from  Equation  5)  accompanying 
increases  in  surface  temperature,  a  destabilizing  mechanism  as  discussed 

further  in  the  next  paragraph. 

In  the  work  described  in  Reference  5,  the  above  equation  for 
as  a  function  of  instantaneous  surface  temperature,  T^,  was  combined  with 
simplified  gas-phase  combustion  model,  a  "full-up"  perturbation  analysis  of 
the  entire  equation  set  describing  surface  and  gas-phase  processes,  and  a 
Fourier  analysis  of  the  subsurface  region  for  development  of  expressions 
relating  the  real  part  of  the  pressure-coupled  response  function 

*  (mj/m)/(P'/P) ]  to  the  dimensionless  angular  frequency  of  oscilla¬ 
tions  (standard  presentation  format)  at  several  values  of  ^ox'^fuel ' 
be  seen  from  Figure  1,  the  differential  dependence  of  oxidizer  and  fuel  abla¬ 
tion  rates  was  predicted  to  have  significant  effects  on  the  real  part  of  the 

f  r  1 

pressure-coupled  response  function,  the  predicted  R  '  ’  peak  increasing 

pc 

dramatically  with  increasing  values  of  that  difference.  Thus,  a  potential 
destabilizing  mechanism  in  AP-composite  propellant  combustion  which  is  not 
present  with  homogeneous  propellants  (and  which  is  not  properly  handled  by  the 
classical  Zeldovich  analysis)  is  suggested.  Most  important,  this  reasoning 
suggests  a  possible  significant  fundamental  approach  to  development  of  a  class 
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Figure  I .  Effect  of  Allowance  for  Dependence  of  Oxidizer/Fuel  Ratio  of  Gases  Leaving  the 
Surface  on  Surface  Temperature.  1091  afsor 


of  AP-composi te  propellants  which  will  burn  stably  over  a  much  wider  range  of 
motor  conditions  than  do  current  formulations;  this  approach  is  centered 
around  development  of  means  of  increasing  the  temperature  sensitivity  of  fuel 
ablation  or  decreasing  that  of  oxidizer  ablation  resulting  in  decrease  (or 

perhaps  even  reversal)  of  the  effect  of  the  gas-phase  heat  release  oscilla¬ 

tions  relative  to  the  surface  temperature  oscillations. 

The  current  program  was  accordingly  directed  at  further  examination 
of  this  approach  to  reduction  of  pressure-coupled  response  functions  of  AP- 
composite  propellants,  to  be  carried  out  in  three  steps: 

1.  Detailed  modeling  involving  perturbation  studies  using  a  complex 
composite  propellant  combustion  model  which  has  been  verified 
against  data,  to  quantitate  potential  effects  of  alteration  of 
oxidizer  or  fuel  ablation  activation  energies. 

2.  Development  of  means  of  tailoring  (changing)  these  ablation 

activation  energies.  At  this  time,  it  appears  that  the  most 
promising  approaches  to  modifying  the  oxidizer  ablation  charac¬ 
teristics  would  involve  occlusion  of  possible  catalysts  in  the 

oxidizer  crystals  or  co-crystallization  of  mixed  oxidizers  (such 
as  ammonium  perchlorate  or  hydroxyl  ammonium  perchlorate,  for 
example).  With  respect  to  modifying  the  temperature  sensitivity 
of  fuel  ablation,  various  approaches  including  the  use  of 
catalysts,  substitution  of  various  moieties  on  the  basic  polymer 
chains,  and  variation  of  crosslinker  type  and  level  are 
available.  In  this  phase  of  the  work,  techniques  for  simul¬ 
taneously  measuring  ablation  rate  and  surface  temperature  of  the 
ingredients  (ablating  alone  under  the  influence  of  an  external 
radiation  flux)  will  be  utilized,  probably  in  collaboration  with 
Prof.  Tom  Brill  of  the  University  of  Delaware. 

3.  T-burner  and  pulsed  motor  testing  with  formulations  containing 
these  modified  ingredients  to  confirm  the  approach. 
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3.0 


STEADY-STATE  MODEL  UTILIZED 


A  complex  steady-state  model  for  combustion  of  composite  solid 
propellants  developed  by  King^^”^^  in  the  late  1970's  and  used  quite  success¬ 
fully  in  predicting  effects  of  various  formulation  parameters  and  pressure  on 
AP/HTPB  composite  propellant  burning  rates  was  selected  as  a  starting  point 
for  this  analysis.  This  model,  based  on  the  same  general  principles  as  the 
classic  Beckstead-Derr-Price  (BOP)  model but  containing  several  major 
modifications,  was  first  further  modified  to  eliminate  a  discrepancy  present 
in  most  (if  not  all)  of  the  BDP-type  models.  [In  these  models,  the  relative 
surface  areas  of  the  oxidizer  and  fuel  are  calculated  via  geometrical  consi¬ 
derations,  with  the  ratio  of  the  surface  area  of  the  fuel  to  the  planar 
oxidizer  surface  area  being  assumed  to  be  equal  to  the  volumetric  ratio  of 
fuel  to  oxidizer  in  the  formulation;  as  a  result,  since  the  volumetric  abla¬ 
tion  rates  of  the  fuel  and  oxidizer  are  in  general  different,  mass  is  not 
conserved  in  most  cases.  That  is,  the  mass  ratio  of  oxidizer  and  fuel  abla¬ 
tion  is  not  in  general  equal  to  the  mass  ratio  of  these  ingredients  in  the 
propellant.]  This  deficiency  has  been  corrected  in  the  modified  steady-state 
model  used  as  a  basis  for  this  study  in  which  an  iterative  approach  regarding 
calculation  of  the  fuel  surface  area  associated  with  a  given  oxidizer  particle 
planar  surface  area  is  utilized  to  ensure  that  the  ratio  of  the  products  of 
areas  and  mass  fluxes  of  oxidizer  and  fuel  _i_s^  equal  to  the  mass  ratio  of  these 
ingredients  in  the  propellant. 

The  basic  model  is  centered  around  an  energy  balance,  the  product  of 
burning  mass  flux  and  energy  requirements  to  raise  ingredients  from  ambient  to 
surface  temperature  (related  to  burning  rate  by  an  Arrhenius  function)  and  to 
vaporize  that  fraction  not  consumed  in  subsurface  reactions  being  equated  to 
the  sum  of  heat  release  rates  from  subsurface  reactions  and  from  two  gas  flame 
zones  (Figure  2).  Thus,  burning  rate  is  controlled  by  three  heat  release 
zones:  (1)  a  thin  subsurface  zone  immediately  adjacent  to  the  surface;  (2)  a 
gas-phase  AP  decomposition  product  monopropellant  flame;  and  (3)  an  extended 
diffusion  flame  in  which  AP  products  and  binder  pyrolysis  products  are  mixed 
and  burned.  Subsurface  heat  release  is  calculated  using  an  estimated  subsur¬ 
face  temperature  profile  substituted  into  a  rate  expression  representing 
subsurface  heat  release  data  measured  by  Waesche  and  Wenograd^^^L  This 
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Figure  2 .  Multiple  Flame  Structure  Model  Schematic  for  APIHTPB  Composite  Propellants. 
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expression  is  integrated  from  the  surface  to  a  depth  where  the  temperature 
equals  the  AP  melting  point  to  obtain  total  subsurface  heat  release.  (This 
procedure  differs  somewhat  from  the  BDP  model,  in  which  subsurface  heat 
release  per  unit  mass  of  propellant  is  assumed  to  be  constant,  independent  of 
burning  rate.) 

For  the  gas-phase,  a  two-flame  approach  was  chosen  (in  contrast  to 
the  three-flame  approach  of  BDP),  the  flames  being  an  AP  monopropellant  flame 
(infinitesimally  thin  flame  sheet  parallel  to  the  propellant  surface)  and  a 
columnar  diffusion  ( Burke-Schumann ^^^^ )  flame.  Three  distances  (FH90,  L^p, 
and  in  Figure  2)  are  important  in  determining  heat  feedback  from  these 
flames.  FH90  is  a  distance  associated  with  90%  mixing  of  the  fuel  and 
oxidizer  gases  from  the  propellant  surface,  while  Lpy  and  L^p  are  reaction 
distances  (products  of  gas  velocity  away  from  the  surface  and  reaction  times) 
associated  with  the  diffusion  and  monopropellant  flames  respectively.  AP 
monopropellant  heat  release  is  assumed  to  occur  at  one  plane,  while  the  diffu¬ 
sion  flame  releases  heat  in  a  distributed  fashion  (as  defined  by  a  Burke- 
Schumann  analysis  described  in  some  detail  in  Reference  6)  between  distances 
and  +  FH90  from  the  propellant  surface.  (Details  of  the  calculational 

procedure  for  obtaining  the  heat  flux  from  these  flames  back  to  the  propellant 

surface  are  also  discussed  in  Reference  6.) 

Details  of  the  equation  development  for  this  model  appear  in 
References  6  and  7.  Included  in  the  model  are  three  "free"  constants  (pre¬ 
exponentials  associated  with  the  subsurface  rate  expression  and  the  two  rate 
expressions  used  to  calculate  the  gas-phase  reaction  times).  Optimized  values 
for  these  constants  are  chosen  using  burning  rate  versus  pressure  data  for 
four  unimodal  oxidizer  AP/HTPB  formulations  and  subsequently  used  foi^  all 

other  calculations  for  noncatalyzed  AP/HTPB  formulations.  Extension  of  the 
model  to  treat  multimodal-oxidizer  formulations  and  metalized  AP/HTPB  formula¬ 
tions  are  discussed  in  References  7  through  9. 

A  brief  description  of  the  steady-state  burning  rate  calculation 
procedure  for  unimodal  oxidizer  propellants  follows.  Again,  the  reader  is 

referred  to  References  6-9  for  more  details  (including  treatment  of  effects  of 
steady-state  crossflow,  which  are  not  considered  in  this  study). 
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1.  Guess  surface  temperature,  T^. 

2.  Calculate  mass  burning  fluxes  and  linear  burning  rates  of 
oxidizer  and  fuel  using; 


ox,s 

(6) 

'fu  ° 

Ap^je'^fu^’^^s 

(7) 

=  m  /p 

(8) 

ox,s 

ox,s  ^ox 

fu  " 

■"fLl'PfU 

(9) 

3.  Utilize  geometrical  arguments  in  combination  with  r^^  to  calcu¬ 
late  the  effective  cured  oxidizer  surface  area  (see 

Reference  7),  ASOX,  of  a  single  oxidizer  particle. 

4.  Calculate  fuel  surface  area  associated  with  a  single  oxidizer 
particle.  In  the  first  loop  of  a  trial-and-error  procedure, 
calculate  this  value  by  equating  the  ratio  of  fuel  surface  area 
to  oxidizer  planar-projection  surface  area  to  the  volumetric 
ratio  of  these  ingredients  in  the  formulation.  In  subsequent 
loops,  calculate  the  ratio  of  [oxidizer  mass  flux  times  oxidizer 
surface  area  divided  by  fuel  mass  flux  times  fuel  area]  to 
overall  oxidizer/fuel  ratio  -  if  this  ratio  is  greater  than 
unity,  raise  the  value  of  fuel  area,  if  less  than  unity,  lower 
the  value.  Continue  this  procedure  until  the  ratio  causes  unity 
and  then  use  a  Newton-Raphson  procedure  to  final  convergence. 

5.  Calculate  the  final  flame  temperature,  T^r,  from  the  input 
oxidizer/fuel  ratio  and  a  table  of  flame  temperature  versus  0/F 
ratio. 

6.  Calculate  energy  transport  and  release  rates  associated  with 
various  processes  in  order  to  perform  an  energy  balance  at  the 
propellant  surface: 

a)  Heat  rate  associated  with  heating  binder  from  the  propellant 
conditioning  temperature  to  the  surface  temperature  at  a 
rate  equal  to  its  ablation  flux. 
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(10) 


b)  Heat  rate  associated  with  ablating  binder  not  consumed  in 
surface/subsurface  reactions 


q,  =  m.  A,  ,Q,  ,  (1-B) 

^3  fu  fuel  fuel, vap  '  ' 


(ii; 


6  =  Fraction  of  fuel  consumed  in  surface/subsurface 
reactions  (discussed  below) 


c)  Heat  rate  associated  with  heating  oxidizer  from  the  propel¬ 
lant  conditioning  temperature  to  the  surface  temperature  at 
a  rate  equal  to  its  ablation  flux 

q,  =  m  A  [C  (T  -T  )  +  Q  1  (12) 

2  ox  SjOx*-  pox^  s  0^  melt, ox-'  ^  ’ 

d)  Heat  rate  associated  with  ablating  oxidizer  not  consumed  in 
surface/subsurface  reactions 


q.  =  m  A  Out (l-a) 
4  ox  s  ,ox  subl '  ' 


(13) 


a  =  Fraction  of  oxidizer  consumed  in  surface/sub¬ 

surface  reactions 


e)  Heat  release  rate  associated  with  surface/subsurface 
reactions 


m  A  OrvnO 
ox  s,ox^EXO 


(14) 


Qexo  “  release  per  gram  of  AP  consumed  in  sur¬ 

face/subsurface  reactions 


The  quantity  a  is  calculated  assuming  that  surface/subsur¬ 
face  reactions  occur  only  in  a  region  where  temperature 
exceeds  the  melting  point  of  AP  and  that  a  reaction  rate 
expression : 
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(15) 


"a  ■ 


is  superimposed  on  an  unperturbed  (unchanged  by  the  heat 
release)  subsurface  temperature  profile  given  by: 


T  -  (T^-T^)  exp  *  Tp 


(16) 


with  a  then  being  calculated  as: 


^  °  ’^AP.melt 

ct  =  f  R  dx/r  (17) 

•'  a  ox,s  '  ' 

X  =  0 


(Details  are  given  in  References  6  and  7).  The  quantity 
6  is  then  calculated  from  a  assuming  a  stoichiometric 
reaction . 

7.  Calculate  the  net  heat  release  per  unit  mass  of  propellant 
associated  with  all  surface/subsurface  processes  as: 


(-q^-q^+q^)  (Wt  Fraction  Oxidizer) 

^*^surf  "  (Oxidizer  MassTlux)  (Oxidi zer ’Surface '^fea)' 

8.  Calculate  rate  of  heat  transfer  back  to  surface  from  the  AP  and 
0/F  flames  depicted  in  Figure  2  using  an  analysis  detailed  in 
References  6  and  7. 

^fdbk  =  '  (LRx.LAp,FH90,m,T,,T^)  (19) 

In  this  analysis,  the  oxidizer/fuel  heat  release  is  assumed  to 
be  a  distributed  release  between  Lp)(  and  +  FH90,  while  the 
AP  flame  heat  release  is  a  sheet  flame  (planar)  heat  release  at 
distance  from  the  surface.  Three  different  calculational 
procedures  are  used  depending  on  whether  L^p  <  +  FH90,  < 

Lap  <  Lpj(  +  FH90,  or  L^p  >  +  FH90.  In  these  calculations, 

Lp)(  and  L^p  are  calculated  as  the  products  of  gas  velocity  off 
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the  propellant  surface  and  reaction  times; 


«x  ■ "  p(n-n^ 

K  V 

I  _  AP  qas,surf  ,, 

'■AP  -  -p(n-ir 

while  tables  generated  from  calculations  performed  using  Burke- 
Schumann  theory  are  used  for  calculation  of  FH90  as  a  function 
of  such  parameters  as  diffusivity,  gas  velocity  off  the  surface, 
flame  stoichiometry  parameteis,  and  oxidizer  and  fuel 
dimensions . 

•  •  •  • 

9.  Compare  the  sum  of  with  plus  if  they  are  not 

equal,  guess  a  new  value  of  surface  temperature  and  return  to 
Step  1.  Repeat  to  convergence. 

10.  After  convergence,  calculate  propellant  mass  flux  and  burning 
rate  from: 


m  A  , 

_ _ ox,s  ox,s _ 1 

propellant  '  (wt  Fraction  Oxidizer)  A  +  A. 

‘  ox,p  fu 

"^propellant  '  "^propel  1  ant ^^propell ant 
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4.0 


COMPARISON  OF  STEADY-STATE  MODEL  PREDICTIONS  WITH  DATA 


As  mentioned  earlier,  there  are  three  "free"  constants,  in 
Equation  16,  and  Kgp  and  K^p  in  Equations  20  and  21,  included  in  the  model 
described  above;  optimized  values  for  these  constants  were  chosen  by  calibra¬ 
tion  against  burning  rate  versus  pressure  data  for  few  unimodal-oxidizer-size 
AP/HTPB  formulations  (three  73-weight-percent  AP  formulations  with  oxidizer 
diameters  of  5,  20,  and  200  microns  and  one  77-weight-percent  AP  formulation 
with  20  micron  diameter  particles).  Burning  rate  data  and  predictions  made 
using  the  optimized  values  of  these  constants  over  a  pressure  range  of  10  to 
150  atmospheres  are  presented  in  Figure  3.  As  may  be  seen,  predictions  and 
data  agree  quite  well  except  for  the  200-micron  AP  formulation,  where 
predicted  values  range  from  5  percent  high  at  10  atmospheres  to  approximately 
40  percent  high  at  150  atmospheres.  Similar  predicted  and  experimental 
burning  rates  versus  pressure  for  five  82-weight-percent  AP  multimodal 
particle  size  formulations  (four  bimodaU  and  one  trimodal),  with  the  same 
values  of  the  constants  being  used  for  the  predictions,  are  presented  in 
Figure  4.  Agreement  between  predictions  and  experimental  data  was  seen  to  be 
excellent  for  all  formulations  except  the  fifty-fifty  20/200-micron  bimodal  AP 
formulation,  where  the  predictions  are  approximately  10%  high  over  the  entire 
pressure  range.  Thus,  it  appears  that  this  steady-state  model,  used  as  a 
basis  for  the  unsteady-state  burning  rate  studies  described  in  the  remainder 
of  this  report,  is  a  good  representation  of  AP/HTPB  composite  propellant 
combustion . 
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Burning  Rate  (cm/sec) 


Figure  3.  Predicted  and  Experimental  Burning  Rate  as  a  Function  of  Pressure  for 

Four  Unimodal  Oxidizer  Formulations  icsi-afsor 
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Burning  Rate  (cm/sec) 


Pressure  (atmospheres) 


Figure  4. 


Predicted  and  Experimental  Burning  Rate  as  a  Function  of  Pressure  for 
Five  Multimodal  Oxidizer  Formulations  ' 
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5.0 


UNSTEADY-STATE  MODELING 


In  this  section,  a  description  of  a  perturbation  analysis  of  the 
steady-state  model  described  in  previous  sections,  consisting  of  a  Fourier 
analysis  of  the  solid  propellant  subsurface  region  plus  a  quasi-steady-state 
analysis  of  the  surface  and  gas-phase  processes  and  including  allowance  for 
perturbations  in  gas-phase  oxidizer  (fuel  ratio),  is  presented.  As  a  critical 
part  of  this  analysis,  the  effects  of  partial  adjustments  of  the  relative 
oxidizer  and  fuel  surface  areas  to  oscillations  in  ablation  fluxes  of  these 
ingredients  were  treated  on  a  parametric  basis.  (For  very  rapid  oscillations, 
there  is  no  time  for  any  surface  area  adjustment,  while  for  slow  oscillations 
approaching  zero  frequency  in  the  limit,  area  adjustments,  required  for  satis¬ 
faction  of  mass  balances  on  oxidizer  and  fuel,  are  total.)  Estimates  of  the 
degree  of  response  of  relative  surface  areas  (between  these  two  limits)  to 
mass  flux  oscillat-ions  were  made  based  on  the  ratio  of  the  oscillation  period 
to  the  time  required  for  burning  of  a  propellant  thickness  equal  to  the 
oxidizer  diameter. 


The  analysis  of  the  response  of  the  subsurface  regions  to  oscilla¬ 
tions  in  heat  feedback  flux  is  a  fairly  standard  one,  used  previously  by  this 
author  (and  many  others)  and  described  in  detail  in  Reference  3.  This 
analysis  results  in  the  following  relationship  between  the  perturbation  values 
of  the  temperature  gradient  just  below  the  surface,  the  surface  temperature, 
and  the  burning  rate; 


where  A 


ai' 

ax 


=  Pp  r  Cp5 


'S 

is  a  complex  function  of  the  oscillation  frequency,  given  by: 


(24) 
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Dimensionless  angular  frequency  of  oscillations 
Solid  propellant  thermal  conductivity 
Solid  propellant  density 
Solid  propellant  specific  heat 

Steady-state  (mean)  burning  rate 

Steady-state  (mean)  burning  mass  flux 
Angular  frequency  of  oscillations  (2nf) 

Axial  distance  (from  surface)  into  propellant 
Perturbation  value  of  temperature 
Perturbation  value  of  burning  rate 


Solution  of  Equation  25,  v;ith  proper  application  of  boundary  conditions  yields 
the  following  expressions  for  the  real  and  imaginary  parts  of  A  as  functions 
of  frequency: 
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real 


+  Wp) 


1/2 


(27a) 


A  =  ^  r  ( 1 

imaginary  2^  ^ ' 
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+  0), 
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(27b) 


Application  of  a  surface  energy  balance  at  any  instant  of  time 
yields  a  relationship  between  the  instantaneous  rate  of  heat  transfer  into  the 
solid  and  the  instantaneous  rate  of  heat  feedback  from  the  gas  phase  above  the 
propellant  surface  (S"^  represents  a  plane  just  above  the  surface,  while  S“ 
represents  a  plane  just  below  it): 


•^g  ax 


ai 
•^s  ax 


m  Q  _ 
s  ^surf 


(28) 


where  Qg^rf  is  the  sum  of  all  heats  of  melting,  vaporization,  and  surface/ 
subsurface  reactions  (sign  convention  chosen  to  be  positive  for  net  exothermic 
process)  per  unit  mass  of  propellant.  (It  should  be  noted  here  that  an  impor¬ 
tant  approximation  has  been  made  in  assignment  of  all  these  heat  release/ 
absorption  processes  to  an  infinitesimally  thin  region  at  the  surface  rather 
than  treatment  of  processes  occurring  over  a  finite  thickness  of  the 


subsurface  region.)  Expansion  of  Equation  28  into  mean  and  perturbation 
quantities  yields: 


1^  il' 
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(29) 


where  the  perturbation  value  of  surface  heat  release  per  unit  mass  of  propel¬ 
lant  is  expressed  as: 


^surf 


dQ 


surf 


dT, 


(30) 


under  an  assumption  that  this  heat  release  is  a  function  of  instantaneous 
surface  temperature  alone. 


Since,  as  discussed  earlier,  a  critical  part  of  the  current  study 
involves  allowance  for  differential  perturbations  in  oxidizer  and  fuel  abla¬ 
tion  rates  with  oscillating  surface  temperature  (due  to  different  activation 
energies  for  the  fuel  and  oxidizer  ablation  processes),  considerable  care  must 
be  taken  in  the  calculation  of  m^  and  as  functions  of  T^.  Two  limiting 

cases  were  analyzed  in  this  study.  In  one  it  is  assumed  that  surface  area 
adjustments  accompanying  mass  flux  oscillations  are  negligible  (high-oscilla¬ 
tion-frequency  limit);  while  in  the  other,  it  is  assumed  that  there  is  time 
for  adjustments  in  fuel  and  oxidizer  surface  areas  to  totally  compensate  for 
the  differential  variations  in  oxidizer  and  fuel  mass  flux  oscillations  such 
that  the  0/F  ratio  of  material  leaving  the  surface  is  constant  at  the  overall 
propellant  0/F  ratio  (low-oscillation-frequency  limit).  (As  will  be  discussed 
later,  parametric  studies  as  regards  frequencies  between  these  two  limiting 
cases  were  carried  out  as  part  of  this  study  via  linear  combinations  of  the 
limiting  case  expressions.)  The  first  limiting  case  (no  surface-area  adjust¬ 
ment)  will  be  referred  to  as  Scenario  A  in  the  remainder  of  the  model  descrip¬ 
tion,  while  the  second  limiting  case  will  be  referred  to  as  Scenario  B.  In 
the  model  developed  in  this  study,  m^  -  T^  and  -  T^  relationships  are 

calculated  numerically  via  calculation  (using  the  steady-state  model  equations 
presented  earlier)  of  m  and  Qj^rf  the  mean  surface  temperature  value  and  at 
slightly  higher  values,  both  with  forbidden  (A)  and  allowed  (B)  surface  area 
adjustments.  With  this  procedure,  values  of  f|y^,  f^g  (where  fj  is  defined 

as  (m'/m)/T'  and  dQ  ,/dT  under  Scenarios  A  and  B  are  calculated. 

'  s  surf  s 
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Combination  of  Equations  28-30  with  consideration 

subsequent  mathematical  manipulation  yields  an  expression  relating  perturba¬ 
tion  values  of  heat  feedback  flux  (from  the  gas  phase)  and  surface 
temperature: 


'’tdbk 
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(31) 


Further  manipulation,  utilizing  a  steady-state  (mean)  energy  balance  at  the 
surface  and  -  fjT^  finally  yields  a  desired  relationship  between  pertur¬ 
bation  mass  burning  flux  and  perturbation  gas-phase  heat  feedback  flux: 
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(32) 


Next,  the  gas-phase  processes  must  be  analyzed  for  development  of  an 

additional  equation  relating  nij.  and  P':  this  equation  is  then 

combined  with  Equation  32  for  elimination  of  yielding  an  expression  for 

the  ratio  of  (mj/ij)  to  (P'/F),  defined  as  the  pressure-coupled  response 
function. 


With  the  gas-phase  analysis  employed  in  this  model  (described  in  the 
steady-state  modeling  section),  there  are  basically  only  three  independent 
parameters  influencing  the  instantaneous  heat  feedback  flux  to  the  propellant 
surface;  these  are  the  instantaneous  surface  temperature  (T^),  the 
instantaneous  flame  temperature  (Tf:),  and  the  instantaneous  pressure  (P).  [Of 
course,  affects  many  other  parameters  which  in  turn  influence  the  heat 
feedback  flux  (e.g.,  fuel  surface  area,  oxidizer  surface  area,  fuel  mass  flux. 
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oxidizer  mass  flux,  fractions  of  oxidizer  and  fuel  consumed  in  surface/ 
subsurface  reactions.)]  As  indicated  earlier,  the  gas  phase  processes  are 
assumed  to  be  quasi-steady-state.  Accordingly,  the  basic  procedure  employed 
in  this  analysis  involves  calculation  of  partial  derivatives  relating 
perturbations  in  heat  feedback  flux  to  perturbations  in  surface  temperature, 
flame  temperature,  and  pressure. 


First,  the  surface  temperature  is  set  equal  to  the  mean  surface 
temperature  value  (calculated  from  the  steady-state  analysis  described 
earlier)  and  the  quasi-steady-state  gas  phase  analysis  (same  as  described 
earlier)  is  used  to  calculate  the  gas-phase  heat  feedback  flux  at  the  mean 
(steady-state)  pressure  and  flame  temperature  values.  Next,  the  surface 
temperature  is  perturbed  slightly  (approximately  1  degree  Kelvin)  upward,  with 
pressure  and  flame  temperature  being  held  at  their  mean  values,  and  new  values 
of  heat  feedback  flux  are  calculated  under  two  scenarios.  (In  Scenario  A,  the 
fuel  and  oxidizer  surface  areas  are  held  at  their  steady-state  values,  while 
under  Scenario  B,  they  are  adjusted  to  cause  the  ratio  of  oxidizer  to  fuel 
gases  entering  the  gas  combustion  zones  to  be  equal  to  the  overall 
oxidizer/fuel  ratio  of  the  propellant.)  From  these  three  calculated  heat 
feedback  fluxes,  the  partial  derivative  of  heat  feedback  flux  with  respect  to 
surface  temperature  at  constant  pressure  and  flame  temperature  is  calculated 
for  each  scenario  from: 
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ax 
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Next,  surface  temperature  and  flame  temperature  are  held  at  their 
mean  (steady-state)  values,  and  pressure  is  perturbed  slightly  upward  for 
calculation  of  a  fourth  heat  feedback  flux  value;  the  partial  derivative  of 
heat  feedback  flux  with  respect  to  pressure  at  constant  surface  temperature 
and  flame  temperature  is  then  calculated  as: 
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AP 


(34) 


22 


Finally,  surface  temperature  and  pressure  are  held  at  their  mean  (steady- 
state)  values,  and  flame  temperature  is  perturbed  slightly  upward  for  calcu¬ 
lation  of  a  fifth  value  of  heat  feedback  flux;  the  partial  derivative  of  heat 
feedback  flux  with  respect  to  flame  temperature  at  constant  pressure  and 
surface  temperature  is  then  calculated  as: 
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aq 


fdbk 
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’fdbk,  T^,  +  AT^,  P  -  '’fdbk,  T^,  T^,  P 


(35) 


Using  linear  superposition  principles  (this  is  a  linearized 
analysis),  perturbations  in  heat  feedback  flux  are  then  related  to  pertur¬ 
bations  in  surface  temperature,  flame  temperature,  and  pressure  by: 


fdbk 


ir  +  +  Np' 


(36) 


At  this  point,  we  have  four  equations  (Equations  30,  32,  36,  and 

m^/ffis-  ^2^s^  perturbation  quantities  (T  5,  T’y  P',  m^, 

and  qf(j|3|()*  One  more  equation  is  required  to  permit  the  desired  calculation 

of  m^  as  a  function  of  P'  (the  true  independent  parameter  in  this  study). 

This  final  closure  equation  is  derived  from  perturbation  analysis  of  the  gas- 
phase  energy  balance,  which  may  be  written  as: 
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(37) 


(38) 


'TOT 


Gas-phase  heat  release  (cal/gm) 

Heat  release  in  surface  processes  (cal/gm) 

Total  available  energy  associated  with  a  "slug"  of 
oxidizer  and  fuel  leaving  the  surface  at  a  given 
instant  of  time  (cal/gm) 
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AHypy  is,  of  course,  dependent  on  the  relative  mass  flow  rates  (mass  flux- 
surface  area  products)  of  oxidizer  and  fuel  leaving  the  surface  at  any  time. 
(ASOX  is  the  instantaneous  surface  area  of  oxidizer,  while  AFU  is  the  corre¬ 
sponding  instantaneous  surface  area  of  fuel.)  A  series  of  calculations  of 
this  available  energy  as  a  function  of  oxidizer/fuel  ratio  over  a  range  of 
interest  covering  practical  AP-composite  propellant  compositions  results  in: 


-2050  +  3600 


"’ox  s^5°^  1 


(39) 


The  quantity  inside  the  brackets,  in  turn,  is  a  function  of  the  surface 
temperature;  the  functionality  obviously  differs  under  the  Scenario  A  (no 
surface  area  adjustment  during  oscillations)  and  Scenario  B  (adjustment  of 
surface  areas  to  maintain  constant  0/F  ratio  gas  "slugs"  leaving  the  propel¬ 
lant  surface)  assumptions.  Again  with  use  of  linearized  analysis  and  defini¬ 
tion  of  the  constant  relating  perturbations  in  the  bracketed  term  to 
perturbations  in  surface  temperature  as: 

f3  »  [  ]'  /!•  (40) 
it  is  easily  shown  that  under  Scenario  B,  f3  »  0,  while  under  Scenario  A, 
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(41) 


where  WFO  is  the  weight  fraction  of  oxidizer  in  the  propellant.  Perturbation 
of  Equation  39,  with  substitution  of  Equation  40  then  yields: 


AH|ot  -  3600  f3T; 


(42) 
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with  subsequent  perturbation  of  Equation  38  then  yielding: 


3600  f3T' 


dQ 


surf 

31“ 


Ts  ■ 


where  f^  »  3600  f^  -  clQ5yj.^/dT^, 


(43) 


Expansion  of  Equation  22  into  steady-state  and  perturbation 
quantities  next  yields  (after  considerable  algebraic  manipulation): 
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which,  in  conjunction  with  Equation  43,  provides  the  final  equation  needed  for 
closure  of  this  analysis.  At  this  point,  we  have  six  equations  (Equations  30, 

32,  36,  43,  44,  and  seven  unknown  perturbation  quantities 

(dfdbk’  ^s’  ^f’  "’s’  ^surf’  ^gas^’  worked  into  the  desired 

expression  for  Rpj,  -  (iti5/mg)/(P' /P)  with  sufficient  effort.  In  the  interest 
of  space,  the  required  manipulation  will  not  be  presented  here;  the  final 
expression  for  the  pressure-coupled  response  function  is: 


(45) 


where  procedures  for  calculation  of  A,  B,  L,  M,  N,  dQ^^^^/dT^,  fj,  and  f^ 
under  either  limiting  scenario  as  regards  surface  topology  adjustment  during 
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oscillations  have  been  previously  discussed,  and  and  are  available 

from  the  steady-state  model.  Given  values  for  all  of  these  parameters,  the 
real  part  of  the  pressure-coupled  response  is  calculated  for  various 
frequencies  using  the  following  procedure. 

First,  the  frequency  is  substituted  into  Equation  26  for  calculation 
of  the  dimensionless  angular  frequency.  Equations  27a  and  27b  are  then  used 
to  calculate  the  real  and  imaginary  parts  of  These  are  then  substituted 
into  Equation  45  for  calculation  of  the  real  and  imaginary  parts  of  RpQ.  From 
calculations  of  various  frequencies,  a  map  of  the  real  part  of  the  pressure- 
coupled  response  versus  frequency  is  generated  for  each  case  examined. 

As  indicated  earlier,  two  limiting  case  scenarios  as  regards  surface 
topology  adjustments  during  pressure  oscillations  were  treated  in  the  equation 
development  for  fj  (appearing  in  the  constant.  A),  f^,  and  L;  in 
the  parameter  studies  described  in  the  remainder  of  this  report,  allowance  was 
made  for  partial  surface  area  adjustments  (intermediate  between  the  two 
limiting  cases  of  zero  adjustment  at  very  high  oscillation  frequency  to  full 
adjustment,  for  maintenance  of  constant  0/F  ratio  of  gas  "slugs"  leaving  the 
propellant  surface,  at  very  low  oscillation  frequency).  A  quantity,  ZFRACT, 
defining  the  fractional  location  of  actual  operation  between  the  two  limiting 
cases  (ZFRACT  -  0  for  Scenario  A,  ZFRACT  «  1.0  for  Scenario  B)  was  utilized  in 
these  calculations,  with  the  values  of  each  of  the  quantities  which  are 
scenario-dependent  being  calculated  as: 

Xj  .  (I-ZFRWT)  .  ZFRACT  ^  (46) 

In  the  first  phase  of  the  numerical  studies,  ZFRACT  was  treated  simply  as  a 
model  input,  with  effects  on  calculated  response  function  versus  frequency 
curves  being  examined  parametrically.  In  the  second  phase  of  the  study,  an 
expression  relating  the  value  of  ZFRACT  to  the  ratio  of  the  characteristic 
oscillation  time  (inversely  proportional  to  the  oscillation  frequency)  and  the 
time  required  to  burn  through  a  thickness  of  propellant  equal  to  one  oxidizer 
diameter  was  utilized.  This  time  ratio,  TIMRAT,  was  calculated  as: 
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TIMRAT 


(47) 


f  p  D 
^prop  ox 

i  -  Propellant  burning  mass  flux  (mean) 

f  -  Oscillation  frequency 

Pprop  “  Propellant  density 

Dqj(  -  Average  oxidizer  particle  diameter 

Included  in  the  user  inputs  to  the  computer  code  in  this  study  phase  were  a 
value  of  TIMRAT  above  which  ZFRACT  -  1.0  (TRATMX)  and  a  value  of  TIMRAT  below 
which  ZFRACT  -  0.0  (TRATMN);  between  these  limits,  it  was  assumed  that  ZFRACT 
could  be  related  to  TIMRAT  by  an  expression  of  the  form: 

ZFRACT  -  Kj  +  K2  In  (TIMRAT)  (48) 

where  Kj  and  K2  were  calculated  using  ZFRACT  -  1.0  at  TIMRAT  -  TRATMX  and 
ZFRACT  -  0.0  at  TIMRAT  -  TRATMN,  leading  to: 


*^2  “  In  (TRATMX)  -  In  (TRATMN) 
Kj  -  -K2  In  (TRATMN) 
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6.0 


RESULTS  OF  PARAHETRIC  STUDIES 


As  indicated  earlier,  there  were  two  major  phases  in  the  numerical 
studies  performed  utilizing  a  computer  code  based  on  the  analyses,  described 
in  the  previous  section,  for  calculation  of  ■'■he  real  part  of  the  pressure- 
coupled  response  functions  of  AP-composite  propellants  versus  oscillation 
frequency.  In  the  first  part  of  the  study,  ZFRACT,  a  parameter  characterizing 
the  fractional  distance  from  zero-surface-topology-adjustment  toward  full 
adjustment  for  maintenance  of  constant  oxidizer/fuel  ratio  in  the  gas  phase 
was  treated  parametrically,  with  results  being  calculated  for  various  user- 
input  values  of  this  parameter  (held  constant  over  the  frequency  range 
studied);  in  this  case,  closed-form  expressions  for  the  real  part  of  the 
pressure-coupled  response  function  in  the  form: 


n{real)  TERM4 

•^pc  “  X  +  A/A  +  Terms 


(51) 


where  TERM4,  TERM5,  and  A  are  constant  over  the  entire  frequency  range  for  a 
given  propellant  and  mean  pressure,  could  be  obtained.  In  the  second  phase  of 
the  study,  where  ZFRACT  is  calculated  as  a  function  of  the  ratio  of  oscilla¬ 
tion  time  to  time  required  for  burning  of  one  oxidizer  diameter  thickness  of 
propellant,  ZFRACT  is  obviously  a  function  of  the  oscillation  frequency;  in 
this  case,  such  a  closed-form  expression  cannot  be  obtained. 


6.1  Variant  I  -  Input  Values  of  ZFRACT 

Calculations  in  this  phase  of  the  study  were  first  made  with  base¬ 
line  values  of  the  oxidizer  and  fuel  ablation  activation  energies  for  defini¬ 
tion  of  the  effects  of  the  ZFRACT  parameter  on  pressure-coupled  response 
functions  of  propellant  in  which  these  activation  energies  have  not  been 
modified  by  chemical  alterations  of  the  ingredients.  Subsequently,  the 
effects  of  different  values  of  these  activation  energies  on  predicted  response 
of  burning  rate  oscillations  to  if;;posed  pressure  oscillations  were  examined. 
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6.1.1 


Calculations  Using  Baseline  Ablation  Activation  Energies 


Typical  plots  of  predicted  values  of  the  real  part  of  the  pressure- 
coupled  response  versus  frequency  are  presented  in  Figures  5-8;  ZFRACT  values 
of  0,  0.2,  0.4,  0.6,  0.8,  and  1.0  were  used  in  these  calculations.  As  may  be 
seen,  the  predicted  magnitudes  of  the  real  parts  of  the  response  function 
decrease  significantly  with  increasing  values  of  ZFRACT;  that  is,  allowance 
for  surface  area  variations  partially  or  totally  compensating  for  the  dif¬ 
ferential  sensitivity  of  oxidizer  and  fuel  ablation  fluxes  to  surface  tempera¬ 
ture  variations  results  in  considerable  stabilization  of  the  burning 
processes.  Another  way  of  viewing  the  result  is  that  it  demoristrates  the 
strong  destabilizing  contribution  of  variations  in  the  oxidizer/fuel  ratio  of 
"slugs"  of  propellant  leaving  the  surface  (accompanying  surface  temeprature 
oscillations)  in  the  absence  of  such  surface  area  variation  compensating 
effects . 


Values  of  TERM4,  TERMS,  and  A  (Equation  51),  along  with  maximum 
values  of  the  real  part  of  the  pressure-coupled  response  function  are  tabu¬ 
lated  in  Table  I  for  82/18  AP/HTPB  formulations  for  various  values  of  ZFRACT, 
pressure,  and  ammonium  perchlorate  particle  size.  (Similar  tables  were 
generated  for  77/23  AP/HTPB  and  73/27  AP/HTPB  formulations,  but  are  not 
presented  here  since  they  add  little  additional  information.)  Again,  it  is 
seen  that  the  maximum  value  of  the  real  part  of  the  pressure-coupled  response 
function  decreases  significantly  as  more  and  more  surface  area  adjustment  is 
allowed  (increased  values  of  ZFRACT).  Careful  examination  of  the  results 
reveals  that  the  maximum  value  of  the  real  part  of  the  pressure-coupled 
response  function,  normalized  by  the  steady-state  burning  rate  exponent,  n, 
correlates  quite  well  with  TERMS  values,  this  normalized  maximum  value 
decreasing  as  TERMS  increases.  (See  Figure  9.)  This  behavior  is  also 
exhibited  for  the  77/23  and  73/27  AP/HTPB  cases.  It  is  of  interest  to  examine 
the  components  of  TERMS  further  to  see  what  factors  lead  to  increased  values 
and  thus  to  decreased  real  part  of  the  pressure-coupled  response  function. 

Comparison  of  Equations  45  and  51  permits  derivation  of  an  expres¬ 
sion  for  TERMS  as  a  function  of  various  parameters: 
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Eqx  -  22000  73/27  AP/HTPB 

Efu- 16900  200jl  AP 

P-10atm 

X3QI  -  .0008  cal/cm  sec  °K 


ZFRACT  -  Fractional  Surface 
Area  Adjustment  (Relative  to  SS) 
During  Pressure  Oscillations 


^0  H 
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Figure  5.  Predicted  Rp^.  (real)  vs.  Frequency  Curves  for  Various  Values  of  ZFRACT. 
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ZFRACT  ••  Fractional  Surface 
Area  Change  (Relative  to  SS) 
During  Pressure  Oscillations 
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Figure  6.  Predicted  Rp^,  (real)  vs.  Frequency  Curves  for  Various  Values  of  ZFRACT. 

31 


(real) 


Frequency  (Hz) 


Figure  7.  Predicted  Rp^  (real)  vs.  Frequency  Curves  for  Various  Values  ofZFRACT. 
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Figure  8.  Predicted  Rp^  (real)  vs.  Frequency  Curves  for  Various  Values  ofZFRACT. 
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Table  I.  TERM4,  TERMS,  A,  and  Maximum  Value  of 

Various  Values  of  AP  Size,  Pressure,  and  ZFRACT. 


82/18  AP/HTPB 

^ox  “ 

22000 

Efy  -  16900 

*^sol  id 

-  0.0008  cal/cm-sec®K 

AP  Diameter 

Pressure 

R(real) 

pc 

(P) 

(atm) 

ZFRACT 

TERM4* 

A* 

TERM5*  Max 

1.0 

10 

1.0 

11.131 

8.3933 

+3.6953 

1.22 

1.0 

10 

0.8 

10.933 

8.2435 

+3.0251 

1.32 

1.0 

10 

0.6 

10.734 

8.0937 

+2.3550 

1.43 

1.0 

10 

0.4 

10.B3B 

7.9439 

+1.6848 

1.62 

1.0 

10 

0.2 

10.337 

7.7941 

+1.0146 

1.84 

1.0 

10 

0 

10.138 

7.6443 

+0.3445 

2.13 

1.0 

40 

1.0 

8.6270 

7.0454 

+2.4818 

1.28 

1.0 

40 

0.8 

8.6679 

7.0789 

+1.9318 

1.37 

1.0 

40 

0.6 

8.7089 

7.1124 

+1.3816 

1.53 

1.0 

40 

0.4 

8.7499 

7.1458 

+0.8316 

1.71 

1.0 

40 

0.2 

8.7908 

7.1793 

+0.2815 

1.90 

1.0 

40 

0 

8.8318 

7.2127 

-0.2686 

2.18 

1.0 

100 

1.0 

7. 6366 

6.7556 

+1.6248 

1.30 

1.0 

100 

0.8 

7.6641 

6.7803 

+1.1354 

1.43 

1.0 

100 

0.6 

7.6916 

6.8050 

+0.6461 

1.58 

1.0 

100 

0.4 

7.6192 

6.8297 

+0.1567 

1.71 

1.0 

100 

0.2 

7.6467 

6.8544 

-0.3326 

1.95 

1.0 

100 

0 

7.6743 

6.8791 

-0.8220 

2.25 

7.0 

10 

1.0 

10.747 

8.3677 

+3.8099 

1.20 

7.0 

10 

0.8 

10.563 

8.2243 

+3.1498 

1.27 

7.0 

10 

0.6 

10.379 

8.0810 

+2.4897 

1.37 

7.0 

10 

0.4 

10.195 

7.9376 

+1.8296 

1.54 

7.0 

10 

0.2 

10.011 

7.7943 

+1.1696 

1.72 

7.0 

10 

0 

9.8265 

7.6510 

+0.5095 

1.99 

7.0 

40 

1.0 

6.8037 

7.1141 

+3.5080 

0.83 

7.0 

40 

7.1479 

+3.0135 

0.88 

7.0 

40 

6.8684 

7.1817 

+2.5191 

0.98 

7.0 

40 

7.2155 

+2.0246 

1.08 

7.0 

40 

7.2493 

+1.5301 

1.19 

7.0 

40 

6.9654 

7.2831 

+1.0356 

1.30 

7.0 

100 

1.0 

4.0828 

6.9305 

+3.4805 

0.53 

7.0 

100 

0.8 

4.1027 

6.9642 

+3.0798 

0.55 

7.0 

100 

0.6 

4.1226 

6.9980 

+2.6792 

0.57 

7.0 

100 

0.4 

4.1424 

7.0317 

+2.2786 

0.59 

7.0 

100 

0.2 

4.1623 

7.0654 

+1.8780 

0.64 

7.0 

100 

0 

4.1821 

7.0991 

+1.4774 

0.73 

-  TERM4/(A  +  A/A  +  TERMS) 
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Table  I.  TERM4,  TERMS,  A,  and  Maximum  Value  of  for 

Various  Values  of  AP  Size,  Pressure,  and  ZFRACT  (Cont'd). 


82/18  AP/HTPB 

Eox  ■ 

22000 

-  16900 

*^sol  id 

.  0.0008 

cal/cm-sec“K 

AP  Diameter 

Pressure 

pc 

(U) 

(atm) 

ZFRACT 

TERM4* 

A* 

TERM5* 

20.0 

10 

1.0 

9.0925 

8.3116 

+4.3698 

0.97 

20.0 

10 

0.8 

8.9549 

8.1858 

+3.7508 

1.01 

20.0 

10 

0.6 

8.8174 

8.0601 

+3.1319 

1.08 

20.0 

10 

0.4 

8.6798 

7.9344 

+2.5129 

1.17 

20.0 

10 

0.2 

8.5423 

7.8086 

+1.8940 

1.30 

20.0 

10 

0 

8.4047 

7.6829 

+1.2750 

1.47 

40 

1.0 

4.1385 

7.2869 

+4.0268 

0.50 

40 

0.8 

4.1573 

7.3201 

+3.6296 

0.52 

40 

0.6 

4.1762 

7.3532 

+3.2323 

0.56 

40 

0.4 

4.1950 

7.3864 

+2.8350 

0.60 

40 

0.2 

4.2137 

7.4195 

+2.4378 

0.63 

40 

0 

4.2326 

7.4527 

+2.0405 

0.66 

Hi 

100 

1.0 

2.8102 

7.1923 

+2.8546 

0.40 

100 

0.8 

2.8221 

7.2228 

+2.5211 

0.42 

100 

0.6 

2.8340 

7.2532 

+2.1875 

0.44 

100 

0.4 

2.8459 

7.2837 

+1.8539 

0.46 

100 

0.2 

2.8578 

7.3142 

+1.5203 

0.49 

100 

0 

2.8697 

7.3446 

+1.1867 

0.52 

Sn 

10 

1.0 

4.2907 

7.8660 

+3.4874 

0.54 

10 

0.8 

4.2889 

7.8627 

+3.0963 

0.56 

!||lH 

10 

0.6 

4.2871 

7.8594 

+2.7051 

0.59 

nnl 

10 

0.4 

4.2854 

7.8561 

+2.3140 

0.63 

10 

0.2 

4.2835 

7.8529 

+1.9229 

0.67 

10 

0 

4.2818 

7.8496 

+1.5317 

0.71 

90.0 

40 

1.0 

3.5945 

8.2683 

+1.4943 

0.59 

90.0 

40 

0.8 

3.5448 

8.1539 

+1.1589 

0.61 

90.0 

40 

0.6 

3.4951 

8.0396 

+0.8240 

0.65 

90.0 

40 

0.4 

3.4453 

7.9252 

+0.4879 

0.69 

90.0 

40 

0.2 

3.3956 

7.8108 

+0.1524 

0.74 

90.0 

40 

0 

3.3459 

7.6964 

-0.1830 

0.79 

100 

1.0 

2.9627 

8.5188 

-0.9098 

0.80 

100 

0.8 

2.8980 

8.3329 

-1.2234 

0.89 

100 

0.6 

2.8334 

8.1470 

-1.5370 

0.98 

100 

0.4 

2.7687 

7.9612 

-1.8506 

1.09 

100 

0.2 

2.7041 

7.7753 

-2.1641 

1.20 

100 

0 

2.6394 

7.5894 

-2.4777 

1.36 

*R(rea1)  _ 
^PC 

TERM4/(A  +  A/> 

+  TERM5) 
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Table  I.  TERM4,  TERMS,  A,  and  Maximum  Value  of 

Various  Values  of  AP  Size,  Pressure,  and  ZFRACT  (Cont'd). 


82/18  AP/HTPB 

i 

>< 

o 

UJ 

22000 

Efy  -  16900 

*^solid 

-  0.0008 

ca1/cm-sec®K 

AP  Diameter 

Pressure 

Max 

pc 

(P) 

(atm) 

ZFRACT 

TERM4* 

A* 

TERMS* 

200.0 

10 

1.0 

3.6739 

7.6769 

+1.7430 

0.60 

200.0 

10 

0.8 

3.7118 

7.6649 

+1.4477 

0.62 

200.0 

10 

0.6 

3.7497 

7.7330 

+1.1623 

0.66 

200.0 

10 

0.4 

3.787B 

7.8111 

+0.8669 

0.69 

200.0 

10 

0.2 

3.8264 

7.8892 

+0.6616 

0.74 

200.0 

10 

0 

3.8633 

7.9673 

+0.2662 

0.81 

200.0 

40 

1.0 

3. 6262 

8.0609 

-0.8644 

0.96 

40 

0.8 

3.6013 

7.9966 

-1.1386 

40 

0.6 

3.4776 

7.9420 

-1.4126 

40 

0.4 

3.4636 

7.8876 

-1.6866 

1.21 

40 

0.2 

3.4298 

7.8331 

-1.9607 

1.36 

40 

0 

3.4069 

7.7786 

-2.2347 

1.63 

200.0 

100 

1.0 

2.9379 

8.4323 

-2.6698 

1.40 

200.0 

100 

0.8 

2.8828 

8.2741 

-2.8299 

1.66 

100 

2.8277 

8.1168 

-3.0899 

1.89 

100 

0.4 

2.7726 

7.9676 

-3.3601 

2.32 

200.0 

100 

2.7174 

7.7993 

-3.6102 

3.09 

200.0 

100 

0 

2.6622 

7.6410 

-3.8703 

4.76 

-  TERM4/(A  +  A/A  +  TERMS) 
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Figure  9. 
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Variation  of  Maximum  Value  of  Real  Part  of  Pressure-Coupled  Response 

Function  (Normalized  by  Burning  Rate  Exponent)  with  TERMS  (Eqn  51). 
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As  may  be  seen  from  this  equation,  decreases  in  dOgu^f/dT^  (more  negative 
values),  f4  (less  positive  values),  and  L  (less  positive  or  more  negative 
values)  all  lead  to  increases  in  TERMS  and  thus  to  decreases  in  the  maximum 
value  of  examination  of  detailed  intermediate  ouputs  of  the  com¬ 
puter  code  indicate  that  such  decreases  in  these  these  parameters  accompany 
increases  in  ZFRACT  (representing  greater  degree  of  surface  area  adjustment 
during  pressure  oscillations).  Changes  in  L  and  dOg^^^/dT^  with  ZFRACT  are 
roughly  equal  in  their  contribution  to  increases  in  TERMS,  while  changes  in  f4 
are  dominant  for  small  AP  particle  size  cases  but  considerably  less  important 

for  large  AP  cases.  Noting  that  f^  -  and  L  - 

constant  flame  temperature  and  pressure),  we  can  see  that  increases  in  ZFRACT 

lead  to  decreases  in  maximum  values  (accompanying  increases  in 
TERMS)  by  causing  decreases  in  the  partial  derivative  of  heat  feedback  flux 
with  respect  to  surface  temperature  (at  constant  T^igj^g  and  pressure)  and 
decreases  in  the  derivative  of  gas-phase  plus  surface  heat  release  with 
respect  to  surface  temperature  (which,  in  turn,  leads  to  decreases  in  the 
derivative  of  flame  temperature  with  respect  to  surface  temperature,  also 
leading  to  lower  increases  in  heat  feedback  flux  accorvpanying  a  given  increase 
in  surface  temperature. 
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6.1.2  Effects  of  Changes  in  Activation  Energies  of  Fuel  and  Oxidizer 
Ablation  Processes 


In  this  phase  of  the  study,  Rf,.  expressions  were  calculated  and 

freall  ^ 

R^c  "  versus  frequency  curves  developed  for  ZFRACT  -  0  (no  surface  area 

adjustment)  and  ZFRACT  -  1  (total  surface  area  adjustment  needed  for  constancy 
of  0/F  ratio  of  gas  “slugs"  leaving  the  propellant  surface)  cases;  results  for 
82/18  AP/HTPB  formulations  with  1,  7,  20,  90,  and  200  micron  diameter  AP  at 
various  pressures  are  presented  and  discussed  here.  In  one  set  of  calcula¬ 
tions  (for  each  pressure  -  AP  size  combination),  the  activation  energy  of  the 
oxidizer  ablation  process  was  held  at  its  baseline  value  of 
22000  calories/mole  with  the  activation  energy  of  the  fuel  being  varied;  in  a 
second  set  of  calculations,  the  fuel  ablation  activation  energy  was  held  at 
its  baseline  value  of  16900  calories/mole,  while  the  oxidizer  ablation  activa¬ 
tion  energy  was  varied.  In  all  cases,  for  either  ingredient,  the  pre-exponen¬ 
tial  in  the  ablation  rate  expression  was  varied  along  with  the  activation 
energy  so  as  to  keep  the  steady-state  burning  rate  at  the  prescribed  pressure 
constant  as  the  activation  energy  varied.  (Obviously,  other  scenarios,  such 
as  holding  the  pre-exponential  constant,  could  have  been  examined,  but  these 
would  have  made  analysis  of  the  results  considerably  more  difficult/ 
confusing.  Further  studies  with  different  treatments  of  the  pre-exponential 

factors  are,  of  course,  possible.)  Results,  in  the  form  of  R^r  expressions 

freal 1  ^ 

and  maximum  R^^  '  values  versus  binder  ablation  activation  energy  at  - 

22000  and  versus  oxidizer  activation  energy  at  E^^gi  -  16900  are  presented  in 

Appendix  A  for  various  AP  sizes  and  pressures  for  ZFRACT  -  0  and  ZFRACT  -  1. 


Plots  of  the  maximum  value  of  the  real  part  of  the  pressure-coupled 
response  are  plotted  against  the  activation  energy  of  the  binder  ablation 
process  at  a  fixed  value  of  oxidizer  ablation  activation  energy  (baseline 
value  of  22000  calories/mole)  for  various  pressure-particle  size  combinations 
for  ZFRACT  -  0  and  ZFRACT  -  1  (the  two  limiting  cases  as  regards  surface  area 
adjustments  during  oscillations)  in  Figures  10-16.  Similar  plots  of  maximum 
Rpc  '  values  versus  oxidizer  ablation  activation  energy  for  the  baseline 
fuel  ablation  activation  energy  of  16900  calories/mole  are  presented  as 
Figures  17-23.  [As  mentioned  earlier,  the  previous  simplified  analysis  of 
References  3-5  indicated  that  the  real  part  of  the  pressure-coupled  response 
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(Figure  10.  Dependence  of  Maximum  Value  of  Real  Part  of  Pressure-Coupled 

Response  on  Binder  Ablation  Activation  Energy  1091  Arson 
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10  n 


Eqx  -  22000 

Efu  -  Variable,  with  pre-exponential  adjusted 
to  give  same  SS  results  at  P  -  40atm 


82/18AP/HTPB 

200p  AP 
P  -  40atm 
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,  -j - , - ^ - , - 1 - 1 - , - 1 
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Figure  1 1 .  Dependence  of  Maximum  Value  of  Real  Part  of  Pressure-Coupled 
Response  on  Binder  Ablation  Activation  Energy 
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Eqx  -  22000 

Efu  -  Variable,  with  pre-exponential  adjusted 
to  give  same  SS  results  at  P  - 1 0atm 
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Figure  12.  Dependence  of  Maximum  Value  of  Real  Part  of  Pressure-Coupled 

Response  on  Binder  Ablation  Activation  Energy  1091  afsor 
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Eox  -  22000 

Efu  -  Variable,  with  pre-exponential  adjusted 
to  give  same  SS  results  at  P  - 1 0Oatm 
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Figure  13.  Dependence  of  Maximum  Value  of  Real  Part  of  Pressure-Coupled 

Response  on  Binder  Ablation  Activation  Energy  losi  afsor 
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14000  16000  18000  20000  22000  24000  26000  28000 
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Figure  14.  Dependence  of  Maximum  Value  of  Real  Part  of  Pressure-Coupled 
Response  on  Binder  Ablation  Activation  Energy 
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1091-AFSOR 


Figure  15.  Dependence  of  Maximum  Value  of  Real  Part  of  Pressure-Coupled 

Response  on  Binder  Ablation  Activation  Energy  iobi-afsoh 
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Figure  16.  Dependence  of  Maximum  Value  of  Real  Part  of  Pressure-Coupled 
Response  on  Binder  Ablation  Activation  Energy 
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1091-AFSOR 


Figure  17.  Dependence  of  Maximum  Value  of  Real  Part  of  Pressure-Coupled 

Response  on  Oxidizer  Ablation  A  ctivation  Energy  1091  afsoh 
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Figure  18.  Dependence  of  Maximum  Value  of  Real  Part  of  Pressure-Coupled 

Response  on  Oxidizer  Ablation  Activation  Energy  ioqi  afsor 
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Efu  - 16900 

Eqx  -  Variable,  wih  pre-exponent al  adjusted 
to  gve  same  SS  results  at  P  -  lOOatm 
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Figure  20.  Dependence  of  Maximum  Value  of  Real  Part  of  Pressure-Coupled 

Response  on  Oxidizer  Ablation  Activation  Energy  iobi-afsor 
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Figure  21.  Dependence  of  Maximum  Value  of  Real  Part  of  Pressure-Coupled 
Response  on  Oxidizer  Ablation  Activation  Energy 
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Figure  22.  Dependence  of  Maximum  Value  of  Real  Part  of  Pressure-Coupled 

Response  on  Oxidizer  Ablation  Activation  Energy  iobi-w^soh 
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Figure  23.  Dependence  of  Maximum  Value  of  Real  Part  of  Pressure-Coupled 

Response  on  Oxidizer  Ablation  Activation  Energy  iobi-afsor 
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decreased  monotonically  with  decreases  in  EqxIDIZER  “  ^BINDER  designated 
as  AE);  this  corresponds  to  monotonic  decreases  of  maximum  with 
i ncreasi ng  fixed  ^gxiDIZER  with  decreasi ng  ^OXIDIZER  fixed 
^BINDER* ^  As  may  be  seen  from  Figures  10-16,  for  ZFRACT  «  0  (no  surface  area 
adjustment,  as  in  the  earlier  simplified  analysis),  the  maximum  values 
decrease  monotonically  with  increasing  EglNDER  (decreasing  jqj2ER"^BINDEr) 
as  in  that  earlier  study.  However,  even  for  the  no-surf ace-topology-adjust¬ 
ment  cases,  the  predicted  effects  of  varying  at  constant  EbiNDER 
(Figures  17-23)  are  not  as  simple;  for  1-micron  and  7-micron  AP  cases,  the 
maximum  ^  values  decrease  with  decreasing  (decreasing 
^OXIDIZER  "  ^binder)  earlier  study,  but  for  20-micron  AP,  these 
values  are  essentially  invariant  with  changes  in  while  for 
90-micron  and  200-micron  AP,  the  maximum  Rp^^^^^  values  increase  strongly  with 
decreasing  EqxidIZER  (decreasing  EqxIDIZER  “  ^BINDER^  direct  contradiction 
to  the  previous  results.  With  complete  surface  area  adjustment  during  oscil¬ 
lations  (ZFRACT  -  1),  giving  constant  0/F  ratio  for  gases  leaving  the  propel¬ 
lant  surface,  the  maximum  value  of  the  real  part  of  the  pressure-coupled 
response  either  increases  or  remains  constant  (depending  on  oxidizer  size  and 
pressure)  with  increasing  binder  ablation  activation  energy  (decreasing 
^OXIDIZER  ■  ^BINDER^  fixed  oxidizer  ablation  activation  energy,  while  it 
decreases  with  increasing  oxidizer  ablation  activation  energy  (increasing 
^OXIDIZER  "  ^BINDER^  constant  fuel  ablation  activation  energy  (opposite  of 
predictions  from  the  old  study).  For  most  cases,  there  seems  to  be  little 
improvement  (decrease)  in  pressure-coupled  response  relative  to  baseline 
values  (marked  on  the  figures)  available  from  modification  of  either  binder  or 
oxidizer  ablation  activation  energies. 

The  maximum  values  of  the  real  part  of  the  pressure-coupled  response 
(maximum  in  the  Rp^^^^^  versus  frequency  curves)  were  found  to  correlate  quite 
well  with  TERMS  in  the  response  function  equation  (Equation  51)  over  the 
entire  range  of  activation  energies  and  ZFRACT  values  examined  for  each  pres¬ 
sure-particle  size  combination  (Figures  24-30),  with  these  maximum  values 
decreasing  with  increasing  values  of  TERMS.  Normalization  of  these  maximum 
values  by  the  burning  rate  exponents  yields  similarly  good  correlation  with 
TERMS;  a  summary  of  these  results  is  presented  in  Figure  31. 
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Figure  24.  Maximum  Value  of  Real  Part  of  Pressure-Coupled  Response  vs.  TERMS 
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Figure  26.  Maximum  Value  of  Real  Part  of  Pressure-Coupled  Response  vs.  TERMS 
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Figure  27.  Maximum  Value  of  Real  Part  of  Pressure-Coupled  Response  vs.  TERMS 
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Figure  28.  Maximum  Value  of  Real  Part  of  Pressure-Coupled  Response  vs.  TERMS 
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Maximum  Value  of  Real  Part  of  Pressure-Coupled  Response  vs.  TERMS 
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Figure  30.  Maximum  Value  of  Real  Pan  of  Pressure-Coupled  Response  vs.  TERMS 
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Figure  31.  Dependence  of  Maximum  Values  of  Real  Part  of  Pressure-Coupled  Response  Divided 
By  BR  Exponent  on  TERMS  for  Several  Pressure-AP  Size  Combinations 
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Identification  of  factors  leading  to  changes  in  TERMS  (and  thus  in 
maximum  value  of  the  real  part  of  the  pressure-coupled  response)  with  changes 

^OXIDIZER  ®  fixed  value  of  EbINDER^  changes  in  EglNDER  ®  fixed 
value  of  EqxIDIZER^  interest.  Examination  of  computer  printout  details 
for  various  cases  in  combination  with  Equation  52  (with  recognition  that 

fj  -  (m'/m)/T^,  L  -  dq^dbk'^^^s  constant  T^.P)  and  f^  -  dOg^^/dT^ 

results  in  such  definition  for  various  pressure-particle  size  combinations, 
both  for  ZFRACT  -  0  and  ZFRACT  -  1,  as  presented  in  the  remaining  pages  of 
this  section.  As  may  be  seen,  several  factors  contribute  to  changes  in  TERM5 
with  changes  in  EqXIDIZER  ^BINDER  other  parameters  being  held 
constant);  the  relative  importance  of  these  factors  is  seen  to  vary  from  case 
to  case,  negating  the  possibility  of  encompassing  generalizations. 
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6.2 


Variant  II  -  Calculations  with  ZFRACT  Being  a  Function  of  a  Charac¬ 
teristic  Oscillation  Time  -  Burn  Time  Ratio 


As  with  the  studies  involving  input  values  of  ZFRACT  (presented  in 
the  last  section),  discussions  of  studies  using  an  expression  relating  ZFRACT 
to  the  ratio  of  characteristic  oscillation  time  to  characteristic  time  for 
passage  of  the  combustion  front  through  a  thickness  of  propellant  equal  to  one 
oxidizer  particle  diameter  are  broken  into  two  parts,  calculations  with  base¬ 
line  values  of  the  oxidizer  and  fuel  ablation  activation  energies,  and 
examination  of  the  effects  of  various  values  of  these  activation  energies  on 
pressure-coupled  response  functions.  In  both  cases,  the  following  relation¬ 
ship  between  ZFRACT  (degree  of  surface  area  adjustment  during  oscillations) 
and  the  time  ratio  (TIMRAT)  was  utilized: 

ZFRACT  -  0.0  for  TIMRAT  <  0.1 

ZFRACT  -  ^  ^  ^  TIMRAT  <10  (53) 

ZFRACT  -  1.0  for  TIMRAT  >  10.0 

(Other  expressions  could  obviously  be  used  in  further  studies,  but  i'"  ■  ''d  not 
offset  the  conclusions  drawn  in  any  major  way.)  As  indicated  earlie’^,  use  of 
an  oscillation-frequency-dependent  expression  for  ZFRACT  of  course  precludes 
development  of  closed-form  expressions  for  the  pressure-coupled  response 
functions  for  any  given  pressure-propellant  combination. 

6.2.1  Calculations  Using  Baseline  Ablation  Activation  Energies 

Plots  of  predicted  values  of  the  real  part  of  the  pressure-coupled 
response  versus  dimensionless  frequency,  wp,  (Equation  26)  are  presented  in 
Figures  32-42  (for  pressures  of  10,  40,  and  100  atm)  for  various  AP/HTPB 
ratios  and  AP  particle  sizes.  Included  on  each  figure  (one  figure  o?r  AP/HTPB 
ratio  -  AP  size  combination)  are  expressions  relating  the  oscillation 
frequency  to  the  dimensionless  value  used  as  the  abscissa  for  each  pressure. 
Since  ZFRACT  decreases  (less  surface  area  adjustment)  with  increasing 
frequency,  these  curves  in  most  cases  have  somewhat  different  shapes  than  the 
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Figure  32.  Real  Part  of  Pressure-Coupled  Response  vs.  Dimensionless  Frequency 
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Figure  35.  Real  Part  of  Pressure-Coupled  Response  vs.  Dimensionless  Frequency 
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Figure  36.  Real  Part  of  Pressure-Coupled  Response  vj.  Dimensionless  Frequency 
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Figure  37.  Real  Part  of  Pressure-Coupled  Response  vs.  Dimensionless  Frequency 
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Figure  39.  Real  Part  of  Pressure-Coupled  Response  vs.  Dimensionless  Frequency 
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Figure  4 1 .  Real  Part  of  Pressure-Coupled  Response  vs.  Dimensionless  Frequency 
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Figure  42.  Real  Part  of  Pressure-Coupled  Response  vs.  Dimensionless  Frequency 
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curves  of  Figures  5-8,  where  ZFRACT  was  held  constant  for  each  curve.  At  the 
lowest  frequencies,  the  response  values  correspond  to  ZFRACT  -  1.0  values  on 
the  previous  graphs;  while  at  the  highest  frequencies,  they  are  closer  to  the 
ZFRACT  «  0  values.  That  is,  these  curves  tend  to  track  across  the  curves 
shown  earlier,  moving  from  the  top  curve  to  the  bottom  curve  as  frequency 
increases  (left  to  right).  (Detailed  computer  printouts  include  ZFRACT  values 
at  each  frequency,  permitting  quantitative  demonstration  of  this  behavior.) 
As  may  be  seen  from  Figures  32-34,  the  pressure-coupled  responses  for  the  1-, 
7-,  and  20-micron  AP,  82/18  AP/HTPB  cases  are  fairly  low,  without  any  sharp 
peaks.  As  the  AP  size  is  increased  to  90  and  200  microns  for  the  82/18  cases, 
the  pressure-coupled  response  versus  dimensionless  frequency  curves  begin  to 
exhibit  significant  peaks  (Figures  35  and  36),  particularly  at  high  pres¬ 
sure.  With  77/23  AP/HTPB  formulations,  increased  sensitivity  to  pressure 
oscillations  is  seen  (Figure  37-39);  this  is  particularly  strong  for  the 
200  micron  AP  case.  In  fact,  for  the  200  micron  AP  case  at  P  =  100  atm,  the 
formulation  is  actually  intrinsically  unstable  (the  dive  in  the  response  to 
very  low  negative  values  is  indicative  of  such  intrinsic  instability). 
Finally,  with  73/27  AP/HTPB  formulations,  the  responses  are  still  higher; 
again,  intrinsic  instability  is  predicted  for  the  200-mi cron-AP,  P  *  100  atm 
case. 

6.2.2  Effects  of  Changes  in  Activation  Energies  of  Fuel  and  Oxidizer 
Ablation  Processes 


Plots  of  the  calculated  real  part  of  the  pressure-coupled  response 
versus  oscillation  frequency  are  presented  in  Figures  43-60  for  various  values 
of  oxidizer  ablation  activation  energy  (at  constant  fuel  ablation  activation 
energy)  and  various  values  of  fuel  ablation  activation  energy  (at  constant 
oxidizer  ablation  activation  energy)  for  several  AP/HTPB  ratio-pressure- 
oxidizer  particle  size  combinations.  The  maximum  values  of  this  real  part  of 
the  pressure-coupled  response,  extracted  from  these  plots,  are  next  plotted 
against  oxidizer  ablation  activation  energy  for  various  AP/HTPB  ratio- 
pressure-AP  size  combinations  in  Figures  61-70,  and  against  fuel  ablation 
activation  energy  in  Figures  71-80.  The  following  discussion  will  be  concen¬ 
trated  on  the  results  depicted  in  Figures  61-80. 
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Figure  53.  Real  Part  of  Pressure-Coupled  Response  vs.  Oscillation  Frequency 
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Figure  54.  Real  Part  of  Pressure-Coupled  Response  vs.  Oscillation  Frequency 
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Figure  56.  Real  Part  of  Pressure-Coupled  Response  vs.  Oscillation  Frequency 
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Figure  57.  Real  Part  of  Pressure-Coupled  Response  vs.  Oscillation  Frequency 
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Figure  58.  Real  Part  of  Pressure-Coupled  Response  vs.  Oscillation  Frequency 
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Figure  61.  Dependence  of  Maximum  Value  of  Pressure-Coupled  Response  on  Oxidizer 
Ablation  Activation  Energy  at  Constant  Fuel  Ablation  Activation  Energy 
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Figure  63.  Dependence  of  Maximum  Value  of  Pressure-Coupled  Response  on  Oxidizer 
Ablation  Activation  Energy  at  Constant  Fuel  Ablation  Activation  Energy 
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Figure  64.  Dependence  of  Maximum  Value  of  Pressure-Coupled  Response  on  Oxidizer 
Ablation  Activation  Energy  at  Constant  Fuel  Ablation  Activation  Energy 
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Figure  65.  Dependence  of  Maximum  Value  of  Pressure-Coupled  Response  on  Oxidizer 
Ablation  Activation  Energy  at  Constant  Fuel  Ablation  Activation  Energy 
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Figure  66.  Dependence  of  Maximum  Value  of  Pressure-Coupled  Response  on  Oxidizer 
Ablation  Activation  Energy  at  Constant  Fuel  Ablation  Activation  Energy 
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Figure  67.  Dependence  of  Maximum  Value  of  Pressure-Coupled  Response  on  Oxidizer 
Ablation  Activation  Energy  at  Constant  Fuel  Ablation  Activation  Energy 
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Figure  68.  Dependence  of  Maximum  Value  of  Pressure-Coupled  Response  on  Oxidizer 
Ablation  Activation  Energy  at  Constant  Fuel  Ablation  Activation  Energy 


109 


1081 -AF  son 


,  n  I  i  I  [  I  I  I 

16000  18000  20000  22000  24000  26000  28000  30000 

Eqx  (cal/molt) 


Figure  69.  Dependence  of  Maximum  Value  of  Pressure-Coupled  Response  on  Oxidizer 
Ablation  Activation  Energy  at  Constant  Fuel  Ablation  Activation  Energy 
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Figure  70.  Dependence  of  Maximum  Value  of  Pressure-Coupled  Response  on  Oxidizer 
Ablation  Activation  Energy  at  Constant  Fuel  Ablation  Activation  Energy 
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Figure  71.  Dependence  of  Maximum  Value  of  Pressure-Coupled  Response  on  Fuel 

Ablation  Activition  Energy  at  Constant  Oxidizer  Ablation  Activation  Energy 
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Figure  72.  Dependence  of  Maximum  Value  of  Pressure-Coupled  Response  on  Fuel 

Ablation  Activation  Energy  at  Constant  Oxidizer  Ablation  Activation  Energy 
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Figure  73.  Dependence  of  Maximum  Value  of  Pressure-Coupled  Response  on  Fuel 

Ablation  Activation  Energy  at  Constant  Oxidizer  Ablation  Activation  Energy 
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Figure  74.  Dependence  of  Maximum  Value  of  Pressure-Coupled  Response  on  Fuel 

Ablation  Activation  Energy  at  Constant  Oxidizer  Ablation  Activation  Energy 
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Figure  75.  Dependence  of  Maximum  Value  of  Pressure-Coupled  Response  on  Fuel 

Ablation  Activation  Energy  at  Constant  Oxidizer  Ablation  Activation  Energy 
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Figure  77.  Dependence  of  Maximum  Value  of  Pressure-Coupled  Response  on  Fuel 

Ablation  Activation  Energy  at  Constant  Oxidizer  Ablation  Activation  Energy 
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Figure  78.  Dependence  of  Maximum  Value  of  Pressure-Coupled  Response  on  Fuel 

Ablation  Activation  Energy  at  Constant  Oxidizer  Ablation  Activation  Energy 
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As  may  be  seen  from  Figures  61-70,  the  maximum  value  of  the  real 
part  of  the  pressure-coupled  response  decreases  with  increasing  value  of  the 
oxidizer  ablation  activation  energy  (at  fixed  fuel  ablation  activation  energy) 
for  virtually  every  case  examined,  the  only  exceptions  being  the  73/27  AP/HTPB 
cases,  which  are  not  really  of  practical  interest  anyway.  This  result  is 
opposite  of  predictions  made  with  the  preliminary  model  (References  3-5) 
developed  in  the  1980  time  period  and  discussed  earlier;  in  that  preliminary 
model  it  should  be  recalled  that  the  pressure-coupled  response  real  parts  were 
predicted  to  increase  with  increasing  oxidizer  ablation  activation  energy  at 
fixed  binder  ablation  activation  energy.  Careful  examination  of  the  computer 
printouts  indicate  that  this  change  comes  about  as  a  result  of  treatment  of 
the  condensed-phase  heat  release  depending  on  various  parameters  (not  treated 
in  the  preliminary  model)  in  the  current  modeling  effort;  these  factors  appear 
to  more  than  compensate  for  the  effects  of  the  relative  ablation  activation 
energies  on  the  degree  of  oscillation  of  oxidizer /fuel  ratio  of  gas  pockets 
leaving  the  propellant  surface,  which  by  themselves  tend  to  lead  to  lower 
response  of  the  propellant  burning  rate  to  pressure  oscillations  with 
decreased  values  of  -  E0j|^|Q£fj.  It  is  also  important  to  note  that 
for  most  cases  (with  the  exception  of  200-micron-AP  cases)  little  benefit  (in 
terms  of  reduced  real  part  of  the  pressure-coupled  response)  is  available  via 
variation  of  the  oxidizer  ablation  activation  energy  at  the  fixed  baseline 
values  of  the  fuel  ablation  activation  energy.  For  the  200-micron-AP  cases, 
considerable  reduction  of  the  maximum  real  part  of  the  response  function  can 
be  obtained  by  increasing  Eq^^jqj^er  (opposite  of  what  was  predicted  with  the 
prel iminary  model ! ) 

Figures  71-80  show  that  the  predicted  maxima  in  the  real  parts  of 
the  pressure-coupled  responses  decrease  with  increasing  values  of  the  binder 
ablation  activation  energy  (at  fixed  oxidizer  ablation  activation  energy)  for 
all  cases  examined,  in  qualitative  agreement  with  the  preliminary  studies  of 
References  3-5.  However,  the  magnitude  of  these  decreases,  particularly  for 
the  more  practical  82/18  AP/HTPB  cases,  are  considerably  less  than  predicted 
in  the  preliminary  studies,  indicating  that  little  benefit  is  likely  to  be 
achieved  through  tailoring  of  the  binder  ablation  activation  energies. 
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7.0 


CONCLUSIONS 


Based  on  the  parametric  studies  with  the  current  model  to  date,  it 
appears  that,  for  most  cases,  little  benefit  in  terms  of  reduced  sensitivity 
of  burning  rate  to  pressure  oscillations  is  likely  to  be  obtained  via 
variation  of  fuel  or  oxidizer  ablation  activation  energies.  However,  it  must 
be  acknowledged  that  the  parametric  studies  to  date  are  fairly  limited  in 
scope;  further  studies  with  simultaneous  variation  of  oxidizer  and  fuel 
ablation  activation  energies  might  reveal  favorable  combinations  for  reduction 
of  pressure-coupled  response  as  might  different  assumptions  regarding  changes 
in  pre-exponential  factors  accompanying  changes  in  actuation  energies. 
(Recall  that  it  was  assumed  in  the  parametric  studies  to  date  that  the 
pre-exponentials  would  vary  with  the  activation  energies  in  such  a  manner  as 
to  give  constant  steady-state  burning  rate  at  the  pressure  of  interest;  other 
scenarios  can  obviously  be  examined  with  this  model,  though  interpretation  of 
the  results  will  be  obviously  more  complex  in  terms  of  what  is  being  held 
constant  in  the  comparisons.)  Thus,  it  may  be  concluded  that  although  preli¬ 
minary  exercise  of  the  unsteady-state  burning-rate  model  developed  in  this 
study  indicates  that  the  proposed  approach  of  reducing  pressure-coupled 
response  functions  via  modification  of  the  kinetics  of  the  oxidizer  and  fuel 
ablation  processes  has  rather  limited  promise,  further  numerical  studies  with 
this  model  might  well  prove  beneficial.  Experimental  studies  should  be 
conducted  to  verify  the  preliminary  conclusions. 
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APPENDIX  A 


Tabulations  of  Rp^.  expressions  and  maximum  values  of  the 
real  part  of  the  pressure-coupled  response  as  a  function 
of  binder  ablation  activation  energy  (at  constant 
oxidizer  ablation  activation  energy)  and  as  a  function 
of  oxidizer  ablation  activation  energy  (at  constant  fuel 
ablation  activation  energy) 


A-1 


82/18  AP/HTPB  lu  AP  P  -  100  atm  ZFRACT  -  0.0 


Eqx  -  22000 


'BINDER 


Rpp  Expression 


Max 


14500 


j  _ 7.51682 

PC  ^  ^  6.73817  _2, 06327 
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j  _ 7.67438 

PC  ^  ^  -0.82536 


19500 


7.84455 

.03181 


A  +  +0.52148 


22000 


8.00872 


A  +  +1.81977 


24500 


8.17270 

Vc  ■  ,  .  Li|548  .3  , 2055 


27000 


8.33655 

7.47233 


A  +  +4.42385 


A-2 


82/18  AP/HTPB  lu  AP  P  -  100  atm  ZFRACT  -  1.0 


Eqx  -  22000 


■BINDER 


RpQ  Expression 


Max 


14500 


7.22121 


A  +  +1.58292 


16900 


, _ 7.53665 

PC  ^  ^  +1.62544 


19500 


7.87856 
.06229  .r 


A  +  ---p^  +1.67042 


22000 


24500 


8.21082 


A  +  +1.71207 


8.54129 


A  +  +1.75447 


27000 


_ 8.87230 

A  +  ^  +1.79739 


A-3 


Eqx  -  16900 

^BINDER  *^PC  Expression  Max 


24000 


22000 


19500 


16900 


15700 


14500 


13200 


12000 


10000 


)  8.27121 

PC  ^  ^  -1.28076 


7.67438 


\  +  -0.82536 


,  _  6.92786 

PC  ^  ^  -0.24807 


6.14555 
^MTTo.  34919 

5.78849 

+0.62634 


5.43114 
k87261  .n 


A  +  ■  -  +0.90209 


„  5.04408 

PC  “  A  .  .1.20543 


4.68673 

?(r555~7 


A  +  +1.49120 


.  _ _  4.09097 

PC  ^  ^  +1.97779 


2.68 


2.26 


1.83 


1.47 


1.33 


1.20 


1.07 


0.96 


0.79 


Eqx  -  16900 

EbindER  '^PC  Expression  Max 


24000 


22000 


19500 


16900 


15700 


14500 


13200 


12000 


10000 


8.02618 

X  4  »2.14615 


7.53665 
6.75566  “, 


+1.62544 


,  6.92779 

A  +  +0.96837 


_ 6.29709 

A  +  +0.28119 


_  5.87809 _ 

X  +  +0.002294 


A  +  ^ 


5.59734 

-0.29833 


_  5.29343 

A  +  -0.64765 


5.01364 
4.49686  n 


X  +  -0.97180 


_ 4.55189 

X  +  -1.51770 


1.25 

1.30 

1.39 

1.52 

1.55 

1.65 

1.79 

1.97 

2.44 


A-5 


82/18  AP/HTPB  7y  AP 


P  -  10  atm 


ZFRACT  -  0.0 


Eqx  -  22000 

^BINDER  *^PC  Expression  Max 


14500 


16900 


19500 


22000 


24500 


27000 


9.62518 


PC  A  t  -1.06186 


)  9.82654 

\  *  +0.50946 


10.04453 


PC  ^  ^  7.82^1  +2.21404 


,  _ 10.25358 

PC  TTTi^nrrj.  84654 


,  10.46331 

A  +  »5. 49289 


j  10.67303 

A  ^  8 •31074  +7^14231 


2.89 


1.99 


1.50 


1.23 


1.04 


0.91 


A-6 


82/18  AP/HTPB  7u  AP  P  -  10  atm  ZFRACT  -  1.0 


Eqx  -  22000 


'BINDER 


Rpp  Expression 


Max 


14500 


,  _ 10.89262 

A  +  »3. 85834 


16900 


,  10.74703 

PC  "  A  .  L3^  ,3.80996 


19500 


,  _ 10.61340 

PC  ^  ^  8.26374  +3  75402 


22000 


10.50816 


A  +  +3.72875 


24500 


A  + 


10.42724 

+3.69603 


27000 


, _ 10.37006 

PC  A  +  +3.66935 


A-7 


82/18  AP/HTPB  7u  AP  P  -  10  atm  ZFRACT  -  0.0 


Eqx  -  16900 


'BINDER 


RpQ  Expression 


Max 


24000 


,  _  10.59092 

PC  ^  ^  8:746U  +0^45475 


22000 


9.82654 


A  +  +0.50946 


19500 


)  8.87220 

PC  ^  ^  §-.?QAi6  +0.57543 


16900 


,  _  7.87791 

X  +  +0.66817 


15700 


14500 


7.41896 


X  +  +0.71451 


6.96012 


A  +  +0.76338 


13200 


,  _ 6.46306 

X  +  +0.81907 


12000 


6.00415 


A  +  +0.87207 


11000 


_  5.62131 _ 

A  +  U^ZiT+o.  91831 


A-8 


82/18  AP/HTPB  7u  AP  P  -  10  atm  ZFRACT  -  1.0 


-BINDER 


Rpp  Expression 


Max 


24000 


11.85748 


A  +  +5.05341 


22000 


_ 10.74703 

PC  ^  §.».36765  +3  80996 


19500 


9.41003 


X  +  +2.24120 


16900 


>  8.02986 

A  +  +0.59830 


15700 


7404999 


A  +  -0.16435 


14500 


j  6.78779 

A  +  -0.92930 


13200 


6.12771 

T.77ro4'T 


X  +  -1.76064 


12000 


j _ 5.52511 

A  +  -2.52811 


11000 


5.03852 


A  +  -3.17608 


-30-32 


A-9 


82/18  AP/HTPB  20u  AP  P  -  10  atm 
Eqx  -  22000 

^binder  he  Expression 


ZFRACT  -  0.0 


Max 


14500 


16900 


19500 


22000 


24500 


27000 


8.23249 


PC  >  4  ?.-.52548  _o. 15744 


8.40476 


A  ^  7.68291  ^1^27501 


)  8.59134 

A  4  MP-li  ^2.83319 


)  8.77074 

A  +  44.33767 


,  8.95089 

A  +  8.18145  ^5^86113 


,  9.13030 

'PC  “  8  34545 

KO  A  +  -  --r-—  +7.37783 


1.94 


1.47 


1.18 


0.99 


0.86 


0.76 


A-10 


82/18  AP/HTPB  20u  AP  P  -  10  atm  ZFRACT  -  1.0 
Eqx  -  22000 

^BINDER  Expression  Max 


14500 


16900 


19500 


22000 


24500 


27000 


9.15933 


PC  “  ^  ^  SATITT 


+4.45799 


j  9.09248 

A  +  +4.36980 


9.04151 


PC  A  +  +4.28309 


9.01372 


PC  X  +  +4.20857 


'PC 


A  + 


9.00756 

-g" 

o • &  0 


+4.13900 


9.02165 


'PC  ,  ,  8^  ,,^03193 


1.01 


1.01 


1.02 


1.03 


1.03 


1.04 


A-11 


Eqx  -  16900 

EbindER  '^PC  Expression  Max 


24000 


22000 


19500 


16900 


14500 


13200 


12000 


11000 


10700 


9.05858 

A  +  +1.52498 


8.40476 


A  +  -  ■  +1.27501 


7.58766 
■6.93582  .n 


+0.96903 


6.73653 
6.15889  .n 


X  +  +0.65713 


5.95291 
r.  44170 


A  +  +0.39053 


5.52827 

)5321 


A  +  +0.23632 


_ 5.13584 

A  +  +0.09731 


4.80997 
:J958i5  „ 


A  +  -0.00182 


4.71206 
■4'.'50622  n 


-0.06620 


1.48 


1.47 


1.46 


1.44 


1.41 


1.40 


1.39 


1.38 


1.39 


A-12 


82/18  AP/HTPB  20u  AP  P  -  10  atm  ZFRACT  -1.0 
Eqx  -  16900 

^BINDER  Expression  Max 


24000 


22000 


19500 


16900 


14500 


13200 


12000 


11000 


10700 


9.99420 


)  9.09248 

A  ^  1.1^1  44.36980 


)  .  7.99147 

A  •!•  42.52256 


6.86146 


PC  X  +  +0.59310 


,  _ 5.86258 _ 

PC  X  +  ~-1.20804 


!  .  5.33201 

A  t  -2.18768 


,  _ 4.84988 

A  t  -3.09135 


, _ 4.44655 

PC  ^  ^  ii.06373  .3^84532 


4.31590 


PC  ^  ^  3.944_18  .4^05345 


0.94 


1.01 


1.17 


1.47 


2.27 


3.65 


11.5 


00 


A-13 


82/18 

Eqx  -  22000 
^BINDER 

14500 

16900 

19500 

22000 

24500 

27000 


AP/HTPB  90y  AP  P  «  100  atm  ZFRACT  -  0.0 


RpQ  Expression 


Max 


2.58714 


PC  ^  .  7.T3368~TTT777 
A  +  — T - 2.83644 


,  _ 2.64197 

-2.48027 


2.70033 


A  +  l-WM  -2.09208 


2.75623 


'PC 


,  _ 2.81130 

7  4  -1.33698 


1.67 


1.37 


1.14 


0.99 


0.88 


'PC 


A  + 


2.86952 

8. '24371 


-0.95540 


0.79 


A-14 


82/18 

Eqx  -  22000 

^BINDER 

14500 

16900 

19500 

22000 

24500 

27000 


AP/HTPB  90u  AP  P  -  100  atm  ZFRACT  -  1.0 


RpQ  Expression 


Max 


•^PC 

^PC 

‘^PC 

“^PC 

'^PC 

*^PC 


3. 03954 

A  +  _o. 48777 

_ 2.96894 

A  +  -0.90409 

2.89115 

A  +  -1.34409 

_ 2.82471 

A  +  -1.76287 

2.76329 

A  +  -2.17554 

2.71645 

A  +  -2.58686 


0.72 

0.80 

0.89 

1.01 

1.19 

1.44 


A-15 


82/18  AP/HTPB  90u  AP  P  -  100  atm  ZFRACT  -  0.0 
Eqx  -  16900 

^BINDER  ^PC  Max 


24000 


22000 


27050 


20000 


19500 


19000 


18750 


,  2.84623 

w  “  rnzipETr:^ 


,  _ 2.64197 

A  +  -2.48027 


,  _ 2.51261 

PC  X  +  ^ -  -3.07004 


,  2.43557 

W  “  ,  ,  Li?895  .3^42534 


,  2.38446 

'pc  ■  ";Op^  .3.65279 


)  2.33304 

A~■^  -3.90127 


,  2.30834 

yc  “  ,  .  6.6Z99I  7~Z^ 
A  +  - r - 4.02083 


0.94 


1.36 


2.00 


2.89 


-4.15 


-7.7 


-15 


A-16 


^BINDER  * 
^BINDER 

24000 

22000 

20750 

20000 

19500 

19000 

18750 


82/18  AP/HTPB  92u  AP  P  -  100  atm  ZFRACT  -  1.0 
16900 

RpQ  Expression  Max 


•^PC 

*^PC 

^PC 

'^PC 

'^PC 


3.29006 


A  +  +0.68064 


2.96894 


A  +  -0.90409 


2.76138 


X  +  -1.89159 


2.64213 


X  +  1^592^  .2.4S854 


2.56038 

A  +  -2.88493 


2.48672 
7.14634 


A  + 


-3.28639 


2.45054 

X  -t-  -3.48661 


0.60 

0.80 

1.06 

1.37 

1.73 

2.44 

3.09 


A-17 


82/18  AP/HTPB  200u  AP 


P  -  10  atm 


ZFRACT  -  0.0 


Eqx  -  22000 

^BINDER  ^PC  Max 


14500 


16900 


19500 


22000 


24500 


27000 


,  3.78394 

k  -  .0.43012 


,  _ 3.86325 

PC  ^  ^  +0.26616 


,  _ 3.94853 

A  +  41.02351 


)  4.03174 

PC  ^  +1.75970 


'PC 


A  + 


4.11418 
8. 484 Jiff 


+2.50210 


'PC 


A  + 


4. 19665 

STESm 


+3.25130 


0.94 


0.81 


0.71 


0.63 


0.58 


0.53 


A-18 


82/18 

Eqx  -  22000 
^BINDER 

14500 

16900 

19500 

22000 

24500 

27000 


AP/HTPB  200u  AP  P  -  10  atm  ZFRACT  «  1.0 


RpQ  Expression 


Max 


^PC 

^PC 

*^PC 

’^PC 

•^PC 

^PC 


_ 3.4S944 

A  ^  LI34_78  ^1^80716 


A  + 


3.67393 


+1.74305 


3.91405 
8.07303 


A  + 


+1.67644 


_ 4.15538 

A  +  61828 


4.40485 

X  \  ^7^^.  ti. 56219 


4.66315 
9.61638 


A  + 


A 


+1.51073 


0.57 


0.60 


0.63 


0.66 


0.70 


0.73 


A-19 


82/18  AP/HTPB  200u  AP 


P  *  10  atm 


ZFRACT  -  0.0 


Eqx  -  16900 

^BINDER  *^PC  Expression  Max 


28000 


25000 


22000 


20000 


18000 


17000 


16800 


,  4.76502 

'PC  “  X  ,  .4.94054 


,  _ 4.31405 

A  4  +2.61423 


,  3.86325 

k  .0.26616 


,  _ 3.56274 

PC  A  4  -1.31178 


,  3.26221 

TC  *  - 6  7277? - 

A  +  -2.89940 


,  _ 3.11168 

A  .  -3.69656 


3.08163 


PC  ^  -3.85623 


0.48 

0.59 

0.81 

1.17 

2.52 

-7.5 

-12 


A-20 


82/18  AP/HTPB  200u  AP  P  -  10  atm  ZFRACT  -  1.0 


Eqx  -  16900 


■BINDER 


RpQ  Expression 


Max  Rpc*''®®') 


28000 


,  _ 4.32354 

X  +  *8.20561 


25000 


3.99156 


X  +  +4.99056 


22000 


3.67393 
’.57685  . 


+1.74304 


20000 


3.46872 

PC  X  +  -0.43942 


18000 


3.27044 


X  +  -2.63742 


17000 


3.17319 


X  +  -3.74268 


16800 


3.15397 


X  +  -3.96397 


A-21 


82/18 

Eqx  -  22000 
^BINDER 

14500 

16900 

19500 

22000 

24500 

27000 


AP/HTPB  200u  AP  P  -  40  atm  ZFRACT  -  0.0 


RpQ  Expression 


Max 


'^PC 

'^PC 

•^PC 

^PC 


3.33682 

X  4  -2.63880 

_ 3.40705 

X  +  -2.23498 

3.48308 

X  +  7.95124  .2  79373 

_ 3.55560 

X  +  SaiT^  -1.36637 

3.62851 

A  +  -0.93601 

3.70110 

A  ^  -0.50280 


1.87 


1.53 


1 

! 

1 

1.28 


1.12 


0.99 


0.90 


A-22 


82/18  AP/HTPB  200u  AP _ P  -  40  atm  ZFRACT  »  1.0 

Eqx  -  22000 

^BINDER  RpC  Expression  Max 


14500 


16900 


19500 


22000 


24500 


27000 


j  _ 3.48070 _ 

PC  7~rZllplI~_o.58492 


,  3  52641 

^pc "  7~7XTj^5~r77rrr 

A  +  - ! - 0.86174 


,  3.59274 

PC  -  TTinorirrrr^ 


,  _  _ 3.65790 _ 

PC  ^  +  8il|0§3  -1.44836 


j  -  3.73296 

A  +  -1.73133 


,  3.81534 

A  >  -2.01144 


0.89 

0.96 

1.06 

1.16 

1.28 

1.42 


A-23 


82/18  AP/HTPB  200y  AP 


P  -  40  atm 


ZFRACT  -  0.0 


Eqx  -  16900 
^BINDER 
28000 

26000 

24000 

22000 

21000 

20000 

19500 

19200 

19000 


Rp0  Expression 


^PC 

'^PC 

’^PC 

^PC 

•^PC 

^PC 

^PC 

^PC 

•^PC 


4.20167 

T*  41.23747 

3.93710 

A  -K  tO. 08942 

_ 3.67207 

A  +  -1.06803 

3.40705 

X  +  -2.23498 

_ 3.27439 

A  ^  -2.82181 

3.14244 _ 

A  ^  -3.41150 

3.07668 

A  -3.70688 

_ 3  03669 

A  »  -3.88448 

3.01082 

A  t  -4.00325 


Max 


0.69 


0.81 


1.04 


1.53 


2.09 


3.46 


5.40 


-8.2 


-12.5 


A-24 


Eqx  -  16900 

^BINDER  Expression  Max 


28000 


26000 


24000 


22000 


21000 


20000 


19500 


19200 


19000 


)  4.46237 

A  +  44.27808 


,  _ 4. 14407 

A  +  +2.57545 


)  3.83277 

A  +  ^Z-SOji  +0.86161 


A  + 


3.52641 

-0.86174 


,  3.38130 

k  ,  7^  .j  — 


3.23571 


A  +  ^  -2.59934 


3.16475 _ 

PC  A  ^  ^-^^^-3. 03551 


3.11649 


A  +  -3.29451 


)  3.09142 

A  ^  -3.47124 


0.48 


0.56 


0.69 


0.96 


1.24 


1.83 


2.45 


3.10 


3.80 


A-25 


Eqx  -  22000 

^BINDER  ’^PC  Exp>^6Ssion  Max 


14500 


16900 


19500 


22000 


24500 


27000 


-9.5 


4.47 

2.97 


2.23 


1.79 


1.50 


82/18  AP/HTPB  200u  AP  P  -  100  atm  ZFRACT  -  1.0 


Eqx  -  22000 


'BINDER 


RpQ  Expression 


Max 


14500 


16900 


2.98627 


\  -2.19964 


2.93577 


A  +  -2.56058 


19500 


,  2.89202 

'pc  -  TTizipzrzJiiTJ 


22000 


,  2.85521 

yC  “  ,  .  8.18648  , 

A  +  - jj - 3.31964 


24500 


2.81901 


A  +  -3.68660 


27000 


?  2.79310 

A  ^  -4.04768 


A-27 


82/18  AP/HTPB  200u  AP  P  -  100  atm  ZFRACT  -  0.0 


Eqx  -  16900 


■BINDER 


Rpf»  Expression 


Max 


28000 


_ 3.28368 

A  +  i^l2418  _i  10756 


26500 


3.12836 

\  ^  Q>.^[841  .1^79077 


25000 


)  .  2.97257 

A  -2.47838 


23500 


j  _ _ 2  81853 

A  +  g-08681  .3  17266 


22000 


,  2.66399 

Vc  -  ,  7^  .3^3,077 


21500 


,  _ 2.61229 

A  ■»•  -4.10533 


A-28 


82/18 

Eqx  -  16900 
^BINDER 

28000 

26500 

25000 

23500 

22000 

21500 


AP/HTPB  200u  AP  P  -  100  atm  ZFRACT  -  1.0 


RpQ  Expression 


Max  RpcI'*®®') 


3.86873 


PC  A  .  .1.89025 


3.63023 


PC  ^  ^  .10>_^876  +0.78685 


TC 


A  + 


3.39111 
9771599 


A 


-0.32583 


3.16733 


PC  ^  -1.43827 


'PC 


A  + 


2.93577 
8.42055 


A 


-2.56058 


'PC 


A  + 


2.86271 
■87211162 


A 


-2.93323 


0.54 


0.62 


0.74 


0.95 


1.42 


1.74 


A-29 


